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EXPLOSIVE LOADING OF METALS 


I. Preface 


The explosive loading of metals, especially as related to shaped-charge 
liners and metals driven by explosives, is a vast field, Thousands of 
documents have been written covering many aspects of the explosive loading of 
materials. However, due to the intense interest in this field by the defense 
establishments of many nations, a great deal of the literature is either 
classified, not available in English, or hidden in obscure government or 
industry documents. Also, certain industries (e.g., drilling, mining, metal 
forming, etc.) engaged in explosive loading, do not readily distribute 
information due to their highly competitive nature. This increases the number 
of proprietary and minimal distribution documents. In short, a great deal of 
information is simply not available to, or cannot be found by, even the more 
diligent researcher. 


I am certainly not aware of all the explosive loading information 
available. Furthermore, I will not even attempt to present all the 
information at my disposal. Instead, I will present basic introductory 
material and cover only the major aspects of the explosive loading of metals 
technology and limit my discussion to shaped-charges and the motion of metals 
driven by explosives. The material and literature covered will be that of 
important historical value and that material found to be most useful to me. 


My sources stem primarily from my association with Ballistic Research 
Laboratory. As a direct consequence of this association, I am probably remiss 
in reporting the results of other government and industrial agencies, 
especially those which do not deal directly with the Ballistic Research 
Laboratory. Also, reference to foreign literature will be minimized. Several 
sources have been omitted as being beyond the scope of this text or merely 
redundant or supplemental to material already included. Nevertheless, a large 
shaped-charge formation and penetration bibliography is given. Again, I do 
not claim that this bibliography is complete. 


Finally, a note of thanks to Lisa Ann Weismiller for the meticulous 
typing of a difficult manuscript. 
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II. The Shaped-Charge Concept 


A cylinder of explosive with a hollow cavity in one end and a detonator 
at the opposite end is known as a hollow charge. The hollow cavity, which may 
assume almost any geometric shape such as a hemisphere, a cone or the like 
eauses the gaseous products formed from the initiation of the explosive at the 
end of the cylinder opposite the hollow cavity to focus the energy of the 
detonation products. The focusing of the detonation products creates an 
intense localized force. This concentrated force when directed against a 
metal plate is capable of creating a deeper cavity than a cylinder of 
explosive without a hollow cavity, even though more explosive is available in 
the latter case. As discussed later, this phenomenon is known in the U.S. as 
the Munroe effect. Further, if the hollow cavity is lined with a thin layer 
of metal, glass, ceramic or the like, the liner forms a jet when the explosive 
charge is detonated. This focusing of the liner material causes a still 
deeper cavity to be formed when the explosive is detonated with the base of 
the liner in contact with a metal plate. The cavity formed becomes deeper yet 
when the explosive charge containing the liner is removed some distance away 
from the plate. This distance, for which an optimum exists, is called the 
standoff distance. Devices of this nature are called lined cavity charges or 
shaped-charges. 


In the following sections, these concepts will be pursued in detail. 
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III, History of Shaped-Charges 


This section will discuss shaped-charge jets or the shaped~charge effect 
resulting from the interaction of an explosive with a metallic liner, 
Further, the discussion will be restricted to conventional shaped-charge 
rounds, e.g., conical, hemispherical and other arcuate metallic liners, and 
Self Forging Fragments (SFF) or Explosively Formed Penetrators (EFP) or P 
charges will only be briefly addressed. Explosive welding and cutting or 
tapping charges will not be considered in any detail. 


The history of shaped-charge development is somewhat ambiguous in that 
the British, Germans and Americans all have made significant claims to the 
early development of lined cavity charges. I understand that currently Dr. 
Manfred Held of MBB, West Germany and H. H. Mohaupt of the U.S. are both 
preparing manuscripts on the history and/or development of Shaped~-Charge 
research. Dr. P, C. Chou [1] and D. F. Kennedy [2] have both drafted 
manuscripts on their version of the history and/or development of 
shaped-charges, The individuals mentioned above are senior experts in the 
field and will certainly provide an excellent account of current and/or 
historical shaped-charge developments. Also, Eather and Griffiths of the 
U.K. [3] have provided a history of the U.K. contribution to the field of 
shaped-charges which includes the contributions of Evans, Ubbelohde, Taylor, 
Tuck, Mott, Hill, Pack and others which I will discuss later. Dr. Held may 
include the discussions of References 1, 2 and 3 in his manuscript. 


Recently, D. F. Kennedy published a description of the shaped-charge 
concept [4] including jet formation, a brief history of the shaped-charge 
concept, jet disturbance mechanisms and modern day applications of 
shaped-charge devices, I highly recommend Kennedy (4] as required reading for 
the novice in the shaped-charge field. 


We will define the term shaped-charge to apply to explosive charges with 
lined or unlined cavities. The cavity is formed in the end of the explosive 
charge opposite the point of initiation. The term shaped-charge, however, has 
@ more general meaning, e.g., Cook [5}. The shaped-charge is also referred to 
as the hollow charge (in the U.K.), the cumulative charge (in the USSR), or 
the Hohlladung (in Germany). 


The hollow cavity effect was first observed by Max von Foerster in 
1883 [6,7] and apparently rediscovered by Charles Munroe in 1888 
[8,9,10,11,12]. The hollow charge or cavity effect is known in the U.S. and 
U.K. as the Munroe effect after Munroe's well documented experiments in 1888. 
Munroe detonated blocks of explosive in contact with steel plates. The 
explosive charge had the initials U.S.N. (United States Navy) inscribed on the 
charge opposite the point of detonation. These initials were reproduced on 
the steel plate. Munroe further observed that when a cavity was formed in a 
block of explosive, opposite the point of initiation, the penetration, or 
depth of the crater produced in the target, increased. In other words, a 
deeper cavity could be produced in a target using a smaller mass of explosive! 
The increase in penetration results from the focusing of the explosive gases 
(detonation products) by the hollow cavity. The depth of the crater in the 
target can be further increased by displacing the hollow charge some optimal 
distance from the target, i.e., increasing the standoff-distance for a lined 
cavity charge. : 
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Early German reference to the hollow cavity effect, excluding von 
Foerster, occurred in 1911 patents in Germany and the U.K. [13,14]. Also, in 
1886, Gustov Bloem [15] of Dusseldorf patented a shell for detonating caps 
(U.S. Patent No. 342, 423) which resembles a shaped-charge with a 
hemispherical liner. 


Early British development of the hollow cavity warhead was reported in 
1913 according to Reference [16]. The earliest known Soviet mention of the 
shaped-charge effect was by Sucharewski in 1925 [17]. The first Italian paper 
on the shapedecharge effect was by Lodati in 1932 [18]. Kennedy [2] provides 
an excellent early history of the shaped-charge effect including evidence that 
this effect may have been known prior to von Foerster using low explosives in 
mining applications such as mudcapping. 


Studies involving the shaped-charge effect continued from 1911 to the 
late 1930's, mostly by German scientists, see, e.g., Kennedy [2]. The 
contributions of Watson [19] on percussion fuzes and Wood [20] on SFF, EFP, 
Misznay-Schardin devices, ballistic dises or P-charge projectiles appeared 
during this period. 


The lined-cavity shaped-charge research accelerated tremendously between 
1935 and 1950 due primarily to World War II and the application of 
shaped-charges to the bazooka, panzerfaust and other devices. The discoverers 
of the lined cavity effect were Franz Rudolf Thomanek for Germany and the 
Swiss, Henry Hans Mohaupt for the U.K. and the U.S. (Unless we credit Bloem 
£15] with this discovery). At any rate, Thomanek and Mohaupt perfected this 
concept and developed the first effective lined cavity shaped-charge 
penetrators. Dr. Thomanek's early work is given in References [21,22,23,24] 
and Dr. Mohaupt's work is given in References [25,26,27,28]. Kennedy provides 
additional references and discussion relating to the work of Mohaupt and 
Thomanek [2]. Kennedy also points out that Munroe used a form of a lined 
eavity charge consisting of an open tin can surrounded by dynamite sticks to 
pierce a steel safe [12]. 


Mohaupt, using lined cavity charges, designed practical military devices 
ranging from rifle grenades, to mortars, to 100 mm diameter artillery 
projectiles. These devices were test fired at the Swiss Army Proving Ground 
at Thun, at Mohaupt's Laboratory and at the French Naval Artillery Proving 
Ground at Gavre. These results were also demonstrated to the U.K. who then 
began development programs of their own. Following the early results of World 
War II, the French Government authorized the release of Mohaupt's information 
to the U.S. and in late 1940 tests were conducted at Aberdeen Proving Ground, 
Maryland using several aspects of lined cavity shaped-charges [28]. The U.S. 
accepted the program, classified it, and thus excluded Dr. Mohaupt from the 
effort but produced the 2.36 in HEAT machine gun grenade and the 75 mm and 105 
mm HEAT artillery projectiles in 1941. Later the machine gun grenade was 
modified to include a rocket motor and a shoulder launcher and became the 
bazooka. The bazooka was first used by the U.K. in North Africa in 1941. 
Other HEAT rounds were fired from tank mounted howitzers [2,28]. 


The military research, introduced by Mohaupt, was supported by private 
firms such as DuPont and the Carnegie Institute of Technology. The Carnegie 
group (C.I.T.) employed some outstanding researchers which contributed much of 
the current shaped-charge knowledge. The leaders at Carnegie were 
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Heine-Geldern, N. Rostoker, Emerson Pugh and his student Robert Eichelberger, 
Eichelberger is the former Director of the Ballistic Research Laboratory, 


Shaped-charge theory continued to develop during the 1950's boosted by 
the Korean War [5,29,30,31,32,33,34,35,36]. During this time period 
tremendous progress was made toward the understanding of the phenomena 
associated with shaped-charge jets. Photographic techniques were employed to 
observe the jet process and analytical models were developed. Efforts were 
made to improve existing shaped-charge liners; to use detonation wave shapers; 
to provide spin compensation via fluted liners; to provide shaped-charge 
follow-through mechanisms; and to enhance the overall system performance, 


Currently, shaped-charge research continues in order to devise successful 
countermeasures to the advanced armors currently fielded and/or contemplated, 
see e.g., Kennedy [4]. Studies which originated in the 1950's still continue; 
notably, torpedo applications of shaped-charge rounds, anti-aircraft rounds, 
fragmentation rounds, multi-staged or tandem liners, long standoff rounds, 
non-conical liners and non-copper liners. Also metallurgical and chemical 
aspects of the liner material was well as methods of liner fabrication remain 
important. 


Starting in the 1950-60's significant shaped-charge developments were 
made possible by the advent of excellent experimental techniques such as high 
Speed photography and flash radiography. Other improvements resulted from the 
transition from TNT to more energetic explosives, 1.e., from TNT to Comp B to 
Octol and then to pressed explosives, notably LX-14. Also, alternate modes of 
initiation (other than point-initiation) and waveshaping techniques have 
provided warhead design improvements. Other advances stemmed from the 
development of large computer codes to simulate the collapse, formation and 
growth of the jet from a shaped-charge liner, These codes provide for the 
most part, excellent descriptions of the formation of the jet. The codes, 
however, are somewhat hindered by the lack of accurate equations of state of 
the liner material and the explosive and by inadequate constitutive 
relationships (i.e., the relationship betweens stress, strain, strain rate, 
and the velocity and velocity gradients). 


The following sections will address the recent developments in more 
detail and will attempt to explain the shaped-charge principles and related 
areas. The information and references given thus far provide a brief 
background of shaped-charge development. Shaped-charge utilization was made 
possible by the discovery of blasting caps (detonators) by Alfred Nobel 
[28,37] in 1867. The explosive reaction initiated by these blasting caps 
could propagate through a column of explosive without the use of confinement. 
This was termed detonation or brisant explosion [28]. The hollow charge 
effect, described earlier, was exploited around 1910 by von Neumann [28] and 
was developed extensively during World War II. 


The Watson patents [19] greatly enhanced the weaponization and further 
applications of the shaped-charge principle. The Watson percussion fuze 
consisted of a parabolic booster charge with a metal lined hemispherical ("or 
arched shield") cavity at the "output" end of the fuze body to intensify the 
effect of the booster charge even over an airgap. This fuze is, in effect, a 
detonator using the shaped-charge principle. 
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Further developments followed by Wood [20] who deseribed what is known 
today as a self forging fragment, explosively formed penetrator, a P-charge or 
a Misznay-Schardin device. 


The shaped=charge principle was clarified and understood as a result of 
the pioneering flash xeray photographs taken by L. B. Seely and J. C. Clark of 
Ballistic Research Laboratory and J. L. Tuck of England, in 1943. Based on an 
analysis of these photographs, analytical models of the collapse of a lined 
conical shaped-charge were developed by Birkhoff [38], Pugh [39], Evans [40] 
and others as indicated in References [1 and 5]. Important contributions were 
carried out by the Carnegie group mentioned previously and continued by the 
Ballistic Research Laboratory (L. Zernow, R. Eichelberger and Associates). 
Other laboratories making important contributions during this time period were 
the Naval Ordnance Laboratory, Maryland (Solem and August), the Naval Ordnance 
Test Station, California (Throner, Weinland, Kennedy, Pearson and Rinehart), 
Picatinny Arsenal, New Jersey (Dunkle), the Standford Research Institute, 
California (Poulter) and others. The liner collapse and jet formation will be 
discussed in detail in a later section. 


A bibliography and account of the weaponization of the shaped-charge and 
similar principles is given by J. Backofen [41-47]. Backofen's bibliography 
is extensive, especially regarding foreign sources. Earlier, WW II results 
and bibliographical information are given in References [48,49,50]. 
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IV. Introduction to Shaped-Charges 


The directional penetration effect observed when a hollow charge is 
detonated in contact with a steel plate is graphically depicted in Figure 1. 
The crater depth is about one-half of the diameter of the hollow of the 
conical cavity. The cavity is produced by high pressure, high velocity gas 
erosion (the Munroe effect). When the hollow cavity is lined with a thin 
hollow metallic or glass cone, the lined charge results in a much deeper 
crater. Furthermore, when the lined cavity charge is displaced from the 
target block some distance (known as the standoff) the penetration increases 
even more, Figure 1 depicts these three cases. 


The increase in penetration resulting from the lined shaped-charge is due 
to the plastic deformation of the liner material under the intense pressures 
and temperatures generated by the detonation of the high explosive. The 
mechanism of jet formation for metallic conical liners with a semi-angle 


(one-half of the conical apex angle) less than 60° is as described below. 


Figure 2 shows a typical shaped-charge warhead like that described above. 
Note that the explosive charge is not cylindrical, but tapered. This removal 
of some of the explosive weight is termed "boattailing" and does not affect 
the jet collapse mechanism. We need only that the detonation wave front be 
Plane and perpendicular to the longitudinal axis-of-symmetry. 


When the detonator is fired, the detonation wave propagates through the 
explosive with the detonation velocity of the particular explosive used. When 
the detonation front reaches the conical liner, the liner is subjected to the 
intense pressures and high temperatures of the front and begins to collapse. 
The collapse is depicted in Figure 3 for the conical lined shaped-charge shown 
in Figure 2. For the position of the detonation wave front shown in Figure 3, 
the upper (apex) region of the cone has collapsed and collided on the 
axis-of-symmetry. This collision results in plastically deformed liner 
Material under tremendous pressure being extruded along the axis-of-symmetry 
(4,40,51]. This extruded material 4s called the jet. When the pressures 
generated exceed the yield strength of the liner material, the liner behaves 
approximately as an inviscid, incompressible fluid. The cone collapses 
progressively from apex to base under point initiation of the high explosive. 
A portion of the liner flows into a compact slug which is the large, massive 
portion at the rear of the jet. 


A rough analogy can be drawn to the effect produced when a sphere is 
dropped into water. During impact and downward motion in the water, the 
cavity formed moves outwards around the sphere and then reverses to collide 
along the axis-of-symmetry. On collision, vertical jets are formed which 
squirt upwards and downwards [40,52]. 


Crude "rules of thumb" were established by Evans [40] in 1950. Namely, 
about 20% of the inner surface of the metal of the cones is used to form the 
jet. The jet diameter is about one twentieth of the diameter of the cone. 
The jet tip velocity is of the order of the detonation velocity of the 
explosive. The jet velocity decay is a nearly liner function down to about 
one quarter of the detonation velocity at the tail (rear) of the jet. The 
velocity of the slug is of the order of one tenth of the jet tip velocity. 
These estimates are only crude order-of-magnitude estimates and are only for 


Iv-1 


moderate apex angle, copper, conical liners. As we will see, there are 
techniques available to provide more accurate estimates of the jet and slug 
parameters. Also note that the jet formation is strongly dependent on the 
liner geometry, liner material, high explosive geometry, confinement geometry 
and material (if confinement is present), the type of high explosive used and 
the mode of initiation. These factors will be addressed later. In any case, 
the goal is to direct and concentrate energy in the axial direction to enhance 
the damage resulting from the hollow charge. 


For any liner design, a proper match between the charge to mass ratio 
(explosive charge mass to liner mass ratio) is critical. If the liner is too 
thick, the energy losses resulting from internal friction and heating of the 
liner walls during the collapse and the energy losses due to spallation of the 
thick liner will reduce the collision velocity below the value necessary for 
jet extrusion. Also, if the liner wall thickness is too thin, directed flow 
is not achieved due to the loss of structural integrity of the liner. If the 
wall is extremely thin, the liner material may undergo vaporization upon 
collision. , 


Shaped-charges with wide angle cones or hemispherical liners show a 
radically different collapse pattern. Hemispherical liners invert (or turn 
inside out) from the pole and as the detonation wave progresses toward the 
base of the liner, the hemispherical liner approximates a conical liner and 
the inverted collapse pattern reverts to that of a cone. In general, the 
hemispherical and large angle conical liners usually result in larger 
diameter, but lower velocity gradient jets. No massive slug, per se, is 
formed. 


Of course, jet effects or controlled fragmentation are not limited to 
conical or hemispherical liners. Extensive use has been made of other liner 
designs (discussed later) as well as linear and circular charges. The liner 
cross sections for linear and circular charges are wedges and semicircular 
configurations, respectively. The jet produced by a linear charge is in the 
form of a thin ribbon and in the case of a torus (semi-circular) charge a 
tubular spray [28]. 


The collapse mechanism of various shaped-charge liners has been verified 
via flash radiography. These novel experimental results lead to the 
formulation of fundamental jet collapse theories. These theorjes will be 
discussed later. 
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Figure 2. Typical Shaped-Charge configuration. 
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Figure 3. Schematic Collapse of a Typical Shaped-Charge. 


V. APPLICATIONS 


Shaped~charges are extremely useful when an intense, localized force is 
required for the purpose of piercing a barrier. The main application is in 
the military arena including torpedos, missiles, high explosive anti-tank 
(HEAT) rounds including hand held (bazooka type) rounds, gun launched (e.g., 
rifle grenades) cannon launched and various bombs. The targets are armors, 
bunkers, concrete or geological fortifications and vehicles, Attacks agains. 
aircraft and spacecraft are possible. Underwater applications (torpedos) are 
possible with the design such that water does not enter the hollow charge 
area. In fact, most warheads of the type described above contain an ogive to 
cover the liner. This ogive acts as an aerodynamic (or hydrodynamic) shield 
while the projectile is in flight, it can provide a housing for impact fuzes 
or guidance, stability and control electronics, and it provides a built-in 
Standoff designed to improve performance, Figure 4 illustrates the concept of 
standoff distance and other nomenclature we will use. The standoff is the 
distance from the base of the charge (or liner) to the target in question, 
Figure 5, from Mohaupt [28] shows an old HEAT artillery projectile. 


The largest shaped-charge application was known as the Beethoven 
Apparatus or Mistel (mistletoe) project [2,53]. The Mistel technique used a 
fighter aircraft mounted piggyback on the top of a large bomber aircraft. The 
unmanned bomber carried a large warhead in its nose. The warhead consisted of 
a 2m (6 foot) diameter, wide angle, conical Shaped-charge, The liner 


probably had a 120° apex angle, was about 30 mm thick, and made of either mild 
steel or aluminum. The warhead weighed 3500 kg (7700 1bs) and the explosive 
alone weighed 1720 kg (3800 lbs). The fighter pilot flew the combination to 
the target, aimed it, released it, then returned to his base. The Germans 
developed this device near the end of WWII and most of them were captured 
intact, 


A shaped-charge with an aluminum liner was used in an attempt to place 
certain identifiable man-made materials into orbit. According to Kennedy [2], 


a 35° included angle aluminum shaped-charge was installed on a multistage 
rocket assembly. The craft was fired into near-space and the warhead was 
detonated in an attempt to project hypervelocity fragments into earth orbit, 
These tests were conducted at Holloman Air Force Base in 1955-56. It was 
never established whether or not any aluminum particles were detected, 


Along the same line, a colleague of mine, S, K. Golaski of the Ballistic 
Research Laboratory and the old Firestone Tire and Rubber Company (now, a part 
of Physics International), developed shaped-charge meteor simulators for NASA 
[54,55]. The objectives of these studies were to obtain the luminous 
efficiency of a meteor like body of known mass, composition and speed during 
re-entry into the earth's atmosphere. The luminous efficiency is the percent 
of kinetic energy of the body which is converted into visible light as 
observed by photographic study of the visible re-entry tail. The study was 
designed to obtain the required mass, composition and speed from pellets 
generated by a specific shaped-charge design flown on a solid state propellant 
vehicle. The vehicle was used to achieve the necessary velocities. Nickel 
and iron were the shaped-charge liner materials. The rocket engines were 
designed to carry the shaped-charge liner above the atmosphere, turn, and 
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allow the detonated shaped-charge to accelerate at hyper-velocity through the 
earth's atmosphere (re-entry). The intent was to simulate a body (meteor) 
re-entering the earths’ atmosphere. A specialized shaped-charge liner design 
was required along with a specialized bi-explosive waveshaping device [55]. 


Many other specialized shaped-charge applications have been pursued by 
the Departments of Defense of several nations. These specialized designs 
included confinement or tamping of the explosive fill, varying the geometry 
and material used for the confinement, varying the geometry or type of 
explosive used, altering the mode of initiation, using explosive lenses or 
more than one type of explosive or an explosive-non-explosive barrier or gap, 
waveshaping or shaping the detonation wave (usually done to insure a uniform 
wave with a short head height, see Figure 4 for definition of head height), or 
varying the standoff distance. Also, significant effects can be achieved by 
varying the liner material (including the use of non-metals such as glass), 
varying the liner thickness, varying the liner diameter, tapering (or causing 
a gradual wall thickness variation either continuously or discontinuously) or - 
varying the liner geometry. The liner geometry variation may utilize the same 
basic geometry, ée.g., varying the conical apex angle, or may employ a 
radically different liner configuration. Other useful liner geometries are 
hemispheres, truncated (from the equator) hemispheres, disc or dish shaped 
(SFF like) devices, tulips, trumpets, dual angle cones, or a combination of 
the above such as hemi-cones or tandem devices. In fact, any arcuate device 
may be used. Also, spin compensated liners may be used, especially when 
associated with spinning warhead applications. Spin compensation (i.e., 
causing the jet to spin enough to compensate for the spin of the warhead) may 
be achieved by metallurgical spin compensation or by the use of fluted liners 
(fluted liners contain raised ridges either on the outside or inside surface 
of the liner to allow the jet to form with a given angular momentum to 
compensate for the rotation of the warhead in flight). Without spin 
compensation, a jet formed from a spinning, detonated warhead will exhibit 
radial inst:bilities (if the rotation speed is high enough) and dispense 
radially, thus reducing the effective penetration capability. Spin and spin 
compensation will be covered in more detail in a later section. 


A unique shaped-charge warhead design was developed by Kennedy and others 
in 1967 to 1970 as reported by Kennedy [2]. This particular missile design 
used a two-stage, tandem liner designed to produce a precursor or pre-jet to 
remove the guidance and control package located in the ogive of the missile. 
Otherwise, the main jet would have to penetrate the seeker package at a short 
(non-optimum) standoff distance. Alignment difficulties in the tandem 
configuration made it necessary to blend the precursor liner, its standoff 
tube and the main conical liner into a single piece. The result was the 
trumpet design configuration which is used today. This design showed improved 
penetration with a relatively small standard deviation in penetration. 


There exists numerous other tales of warhead development, the point being 
that most of the warhead concepts being pursued today are not original 
concepts. Other warhead/liner designs will be discussed later. 


Another application of shaped-charges is in demolition work. This area 


has both military and industrial application. Buildings, bridges and other 
structures are the common demolition targets. Usually demolition charges are 
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hand placed. Demolition charges are sometimes constructed and placed in sites 
using a plastic explosive and a collection of liners. This technique allows 
for more flexibility and adjustment to the conditions at hand than a 
collection of fixed charges. Of course, experience on the part of the blaster 
is required and this technique (the adjustment of the charge to the task) has 
delayed -bulk charge production and research towards this particular 
application. The shaped-charge principle is also used in construction work to 
break, crack or drill holes in rock. A technique known as mudeapping is 
sometimes used to break rock and usually utilizes an unlined hollow charge. 
Shaped-charges are also used in construction as earth movers, in tunneling, or 
to assist in well drilling. ‘ 


Hollow charges are also used as a source of earth waves for geophysical 
prospecting and seismic exploration, Other applications occur in mining 
(surface or underground), submarine blasting, timber cutting, breaking log 
jams, breaking ice jams and steel mill furnace tapping. 


A jet tapper is a shaped-charge with a conical metal or glass liner. The 
steel mill furnance tapping problem is depicted in Figure 6 [37]. The tapping 
problem requires a means of starting the flow of molten steel once the tap 
hole has been plugged. Steel mill jet tapping is sometimes called salamander 
blasting since salamander is the term used to describe the large mass of iron 
deposited at the base of the blast furnance after it has been in operation for 
some time. The tap hole is dug out as far as possible and jet tapper is 
inserted into the tap hole using a loading pole. The jet tapper is then 
detonated, clearing the plug, usually with a high rate of success (371. 

A commercial, DuPont type jet tapper is schematically illustrated in Figure 7. 


The shaped-charge assembly may be aligned on axis or offset up to 10°, as 
shown in Figure 7, to reduce the line-of-sight thickness of the plug that must 
be perforated. , 


Another industrial application of the shaped-charge is in the internally 
coned end of certain detonators. This indented, lined cavity acts to 
concentrate the effect along the axis, recall Bloem [15]. 


Cook [56] used the original Munroe effect to engrave or stencil letters 
onto metal plates. 


Other examples are in the oil well industry where large diameter, but 
extremely short, lined shaped-charges are used to penetrate various geological 
formations to increase the flow of oil. Oi1 well perforation problems present 
extremely difficult design problems due to the minimal amount of allowable 
space available in the well, the short standoff distances required and the 
hostile environment within the well [57]. 


Linear (wedge shaped, V-shaped or W-shaped) shaped-charges are also used 
as cutting charges, They generate a ribbon shaped jet used to cut metals and 
other materials. Commercial cutting charges are available from such places as 
Ensign-Bickford, Explosive Technology, DuPont and others. For cutting charge 
applications, it is sometimes advantageous to use homemade cutting charges 
which can be optimized to the particular problem on hand. Cutting charges and 
hollow charges are also used as explosive separation devices, as bolt cutters 
and for other applications. The shaped-charge effect is used on systems for 
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separation, deployment, safety destruct and similar spacecraft components. 
Gareia [58], Brown [59] and Sewell [60] provide additional data on the linear 
cutting charge concept. Leidel [61] presents a discussion of the design of a 
annular cutting charge. 


An interesting paper by Jones [62], discusses the perforation of artic 
sea-ice by large, demolition type, shaped-charges. The following table was 
extracted from Reference (62] and illustrates the charge sizes used as well as 
the penetration depth achieved in various targets. 


General references on application are found in Kennedy [2], Cook [5], The 
Blaster's Handbook [37], and Bawn and Rotter [63]. Certain specific 
applications of interest are found in Cook [56], on engraving; Torrey [64] on 
a review of shaped-charge weapons of WW II; Hutt] [65] on underground 
blasting; Clark [66] on concrete fragmentation; Lawrence [67], and Austin, et 
al, [(68,69,70) on boulder breaking, drilling and seismic exploration; Draper, 
et al, (71] on drilling; and Davidson and Westwater [72] on cratering and 
boulder breaking. Excellent overviews are available in Murphy [73] and Baum, 
et al. [74]. Extensive studies are currently underway in cratering (earth 
removal), boulder breaking, and penetration of concrete and geological 
formations at several installations, notably the Waterways Experimental 
Station, Vicksburg, Mississippi. The demolition and cratering work performed 
at the Waterways Experimental Station is too extensive to cover here, but 
their studies cover many of the applications cited earlier. See Joachim [75] 
as one example. 


Additional applications relating to explosive metal interactions, and 
which require an understanding of the jet formation phenomena, are explosion 
welding, explosion calding, or explosion forming of metal parts. The 
principles of explosive welding are beyond the scope of this text. Basically, 
the technique involves using explosives to form bimetallic or trimetallic 
parts. A tremendous savings in material usage is thus possible because a thin 
layer, of an often costly material, can be applied to a cheaper, load bearing 
substrate. The thin noble and expensive material is usually used to provide 
an inert or corrosion resistance barrier to the substrate material. 
Applications occur in chemical vessels, heat exchanges, desalination plant 
piping and other areas where caustic environments exist. With explosive 
welding techniques even metallurgically incompatible metals and immiscible 
‘ metal combinations can be bonded, as well as compatible metal systems. 
References [76, 77 and 78] provide an excellent introduction to the various 
aspects of explosion welding. Later a brief discussion of this important 
field will be provided. 
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Figure 7. The Jet Tapper Charge (Offset) [37] . 
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TABLE I. 
FROM J. M JONES [62] 


DESIGN AND PERFORMANCE DETAILS OF STANDARD DEMOLITION SHAPED CHARGES. 


No. 1, 6 in, No. 11, 30 lb. M2A3 M3 
Total weight (1b.) approx, 10 50 15 40 
Explosive weight (lb.) 6.7 27.5 11.5 30 
Charge diameter (in.) 6.0 12 6.0 9.5 
Charge length (in.) 7.5 11.75 9.44 15.5 
Design Standoff (in.) 5.7 19.0 5.5 15.0 
Casing material steel steel fibre steel 
Cone diameter (in.) 6.0 12.0 4.88 9.0 
Cone angle (deg.) 80 60 60 60 
Cone material steel steel glass steel 
Cone thickness (in.) «140 2225 «350 .150 
Cone weight (1b.) - - 1.7 5.6 
Penetration 
Armour depth (in.) 6 - 12 20 + 
Average diameter (in.) - - 1.5 2.5 
Mild steel depth (in.) 9 - - - 
Concrete depth (in.) 30 72 30 60 
‘Entrance diameter (in.) - 6-8 3.5 

Minimum diameter (in.) 2 - 2 2 
Fresh water ice depth (in.) - , - 72 144 

diameter (in.) - - 3.5 6 
(standoff 42 in.) 
Permafrost depth (in.) - - 72 72 


(30 in. (50 in, 
standoff) standoff) 


diameter (in.) - - 6 - 1.5 8-5 
Crater diameter at - - 28 28 


design standoff (in.) ye19 


11. 


12. 


13. 


REFERENCES 
(Sections I ~ Vv) 


P. C, Chou and W. d. Flis, "Recent Developments in Shaped Charge 
Technology," Private correspondence, April 1984 (Manuscript prepared for 
M. Held, MBB, Germany, and presented at MBB Schrobenhausen, West Germany, 
September 1983). 


D. R. Kennedy, "The History of the Shaped Charge Effect, The First 100 - 
Years," presented at MBB Schrobenhausen, West Germany, September 1983. 


R. Eather and Neill Griffiths, "A U. K. Note on the History of Shaped 
Charges," Royal Armament Research and Development Establishment Report, 
August 1983, presented September 1983 at MBB, Germany. 


Donald R. Kennedy, "The Infantryman VS the MBT," National Defense, ADPA, 
Pp. 27-34, March 1985. 


Melvin A. Cook, The Science of High Explosives, American Chemical Society 
Monograph Series, Reinhold Publishing Corp., NY, 1958. 


Max von Foerster, "Versuche mit Komprimirter Schiessbaumwolle," 1883, 
Mittler and Son, Berlin. 


Max von Foerster, "Experiments with Compressed Gun Cotton," Van 
Nostrand's Engineering Magazine, Vol. 31, July - December 1855, Pp. 
Ti3-119. 

Charles E. Munroe, "On Certain Phenomena Produced by the Detonation of 


Gun Cotton," Newport Natural History Society, Proceedings 1883-1888, 
Report No. 6, CER 


Charles E. Munroe, "Wave-like Effects Produced by the Detonation of Gun 
Cotton," American Journal of Science, 36, 1888, pp. 48-50. 


Charles E. Munroe, "Modern Explosives,'" Seribner's Magazine, Vol. 111, 
NY, January - June 1888, pp. 563-576. 


C. E. Munroe, Executive Document No. 20, 53rd Congress, tst Session, 
Washington, DC, 1894, . 


Charles E. Munroe, "The Applications of Explosives," Popular Science 
Monthly, in 2 parts pp. 300-455, Vol. 56, 1900. 


"Verfahren zur Herstellung von Sprengkorpern," Westfolisch - Anhaltische 


358, 1910/1911, 


"Improvement in Explosive Charges or Bodies," Patent No. 28,030 (UK) to 
Westfolisch - Anhaltische Sprengstoff Actien Gesellschsft (WASAG), 
Berlin, 12 December 1911. 


V-11 


16. 


17. 


20. 


21. 


22. 


23. 


24, 


25. 


26. 


27. 


28. 


29. 


30. 


Gustav Bloem of Dusseldorf, Prussia, Germany, "Shell for Detonating 
Caps," U.S. patent 342, 423 dated May 25, 1886. 


LT H. W. Kline, "The Cavity Charge, Its Theory and Applications to the 
Opening of Explosive Filled Ordnance, etc.," U.S.N.R., Ordnance 
Investigation Laboratory, Naval Powder Factory, Indian Head, Maryland, 15 
August 1945. 


Sucharewski, Technica j Snabschenie Krasnoi Armii, 1925, No. 170, pp. 
13-18, also Woina i Technica, 192, No. 253, pp. 18=24,. 


C. Lodati, "An Explanation of the Explosive Behavior of Hollow Blocks of 
Compressed TNT," Giornale di Chim. ind. ed. appl., Vol. 14, 1932, pp. 130. 


Charles P. Watson for "Percussion Fuzes," 1925 U.S. Patents 1, 524, 011 


and 1, 534, 012 filed 22 September 1921 and 27 August 1923, both issued 
on 14 April 1925. 


R. W. Wood, "Optical and Physical Effects of High Explosives," 
Proceedings of the Royal Society (London), Vol. 157A, 1936, pp. 249-261. 


Heinz Freiwald, "The History of Hollow Charge Effect of High Explosive 
Charges," with forward by Hubert Schardin, for German Academy of Aviation 
Research, 15 September 1941. 


"High Explosive Charge," Hungarian Patent 134378 assigned to Engineers 
Carl Brandmayer and Franz Rudolf Thomanek, Berlin, December 9, 1943. 


Hubert Schardin, "Development of the Shaped Charge," published in 
Wehrtechnische Hefte, Heft 4, 1954. 


Franz Rudolf Thomanek, "Meine Hohlladungs - Aktivitaten," for Periods 
1932-1935, 1938-1945, 195721969 and from 1975 to date, 10 March 1978. 


"An Improved Explosive Projectile," Patent, Commonwealth of Australia, 
assigned to Berthold Mohaupt, Henry Mohaupt and Erick Kauders of France, 
August 1941 (date claimed for patent 9 November 1939). 


"Projectile," U.S. Patent No. 2,419,414 assigned to Henry H. Mohaupt, 
April 22, 1947 originally filed October 3, 1941. 


"Projectile," U.S. Patent No. 2,974,595 assigned to Henry H. Mohaupt, 
1961, filed September 11, 1947. 


Aerospace Ordnance Handbook, edited by F. Pollard and J. Arnold, Chapt. 
11, "Shaped Charges and Warheads," H. Mohaupt, Prentice-Hall, Inc., NJ, 
1966. 


H. Kolsky, C. I. Snow and A. C. Shearman, "A Study of the Mechanism of 
Munroe Charges, Part I - Charges with Conical Liners," Research 
Supplement, 2-2, pp. 89-95, London, 1949. 


H. Kolsky, "A Study of the Mechanism of Munroe Charges, Part II = Charges 
with Hemispherical Liners," Research Supplement, 2-2, pp. 96-98, London, 
1949. 


31. 


32. 


33. 


34. 


35. 


36. 


37. 


38. 


39. 


40. 


44, 


42, 


43. 


yy, 


45. 


46. 


W. M. Evans and A. R. Ubbelohde, "Formation of Munroe Jets and Their 
Action on Massive Targets," Research Supplement, 3-7, London, 1950. 


W. M. Evans and A. R. Ubbelohde, "Some Kinematic Properties of Munroe 
Jets," Research Supplement, 3-8, London, May 1950. 


E. M. Pugh, R. V. Heine-Geldern, S. Foner and E. C. Mutschler, "Kerr Cell 
Photography of High Speed Phenomena," J. of Appl. Phys., Vol. 22, No. 4, 
April 1951, pp. 487-493, 


W. S. Koski, F. A. Lucey, R. G. Shreffler and F. J. Willig, "Fast Jets 
From Collapsing Cylinders," J. of Appi. Phys., Vol. 23, No. 12, December 
1952, pp. 1300-1305. 


Franz R. Thomanek, "Die Erste Hohlladungswaffe - The First Hollow Charge 
Weapon," Explosivstoffe, No. 1, 1959, pp. 9=11. 


Franz R. Thomanek, "Die Entwicklung der Ausgekleiden Hohliadung - The 
Development of the Lined Hollow Charge," Explosivstoffe, No. 8, 1960, pp. 
1T7+179. 


The Blaster's Handbook, E. I. du Pont de Nemours and Co., Wilmington, DE, 
16th Edition, . 


G. Birkhoff, D. MacDougall, E. Pugh, G. Taylor, "Explosive with Lined 
Cavities," J. of Appl. Physics, Vol. 19, No. 6, June 1948. 


E. Pugh, R. Eichelberger and N. Rostoker, "Theory of Jet Formation by 
Charges with Lined Conical Cavities," J. of Appl. Physics, Vol. 23, No. 
5, May 1952. 


W. M. Evans, "The Hollow Charge Effect," Bulletin of the Institution of 
Mining and Metallurgy, No. 520, March 1950. 


Joseph E. Backofen, "The Weaponization of Shaped Charge Technology," 
Proceedings of the 4th International Symposium on Ballistics, 17-19 


October 1978. 


J. E. Backofen, "Shaped Charges Versus Armor," Armor, July-August 1980. 


J. E. Backofen, "Shaped Charges Versus Armor - Part II," Armor, September 
- October 1980. 


J. E. Backofen, "Shaped Charges Versus Armor - Part III," Armor, November 
- December 1980. 


J. E. Backofen and L. W. Williams, "Antitank Mines," Armor, July - August 
1981. 


J. E. Backofen and L. W. Williams, "Antitank Mines - Part II," Armor, 
September - October 1981. 


V-13 


47. 


48. 


49. 


50. 


51. 


52. 


53. 


54. 


55. 


56. 


5T. 


58. 


59. 


60. 


61. 


62. 


J. E. Backofen and L. W. Williams, "Antitank Mines - Part III," Armor, 
November ~ December 1981. 


“Annotated Bibliography of NDRC Technical Reports and Memorandums of 
Division 2, including pertinent division 8 Reports," May 1, 1945, NRDC 
(National Defense Research Council) Memorandum No. A-106M, OSRD Report 
No. 4830B. 


Hill, Mott and Pack, A.R.D. Theoretical Research Report No. 2/44 (January 
1944) and 13/44 (March 1944), 


*protection Against Shaped Charges," Carnegie Institute of Technology, 
NDRC Report No. A-384 and OSRD Report No. 6384, February 1946. 


Marvin E. Backman, "Terminal Ballistics," Naval Weapons Center Technical 
Publication 5780, China Lake, CA, February 1976. 


A. M. Worthington, A Study of Splashes, Longmans Green, London, 1908. 


R. Coles and P. L. Rickson, "Mistletoe - The Deadly Parasite," Air 
Classics Quarterly Review, Vol. 4, No. 3, pp. 38-47, 86-99, Fall 1977. 


"Development and Testing of Advanced Shaped Charge Meteoritic Simulators 
Part II - Calibration of Flight Guns," R. L. Woodall and E. L. Clark, 
Firestone Tire and Rubber Company, NASA CR-66216, November 1966. 


"Design, Development and Testing of Bi-Explosive Shaped Charge Meteoritic 
Simulators," The Firestone Tire and Rubber Company, NASA CR-66615, April 
1968. 


M, Cook, Research, 1948, 1, pp. 474. 


J. S. Reinhart and R. D. Cocanower, "Concerning the Design of an 
Effective Shaped Charge for Oil Well Perforating," J. Appl. Phys., Vol. 
30, pp. 680— 682, May 1959. 


M. A. Garcia, "The End-Initiated, Linear Shaped Charge: An Analytical 
Model," Naval Missile Center, TM-67-64, December 1967. 


G. E. Brown, "A New Optimization Theory for the End-Initiated Linear 
Shaped Charge," M. S. Thesis, Naval Postgraduate School, October 1969. 


Robert G. S. Sewell, "Effects of Velocity and Material Properties on 
Design Limits for Linear Shaped Charges," NOTS Technical Publication 3894 
and NAVWEPS Report 8793, August 1965. 


David J. Leidel, "A Design Study of an Annular - Jet Charge for Explosive 
Cutting," Doctoral Dissertation, Drexel University, Philadelphia, PA, 
June 1978. 


J. M. Jones, "The Design of a Large Shaped Charge Suitable for the 


Perforation of Artic Sea-Ice," Defense Research Establishment, Valcartier 
Report, DREV TN-1907/71, PROJ D97-01-19, May 1971, CAD 913074). 


yv=14 


63. 


64. 


65. 


66. 


67. 


68. 


69. 


70. 


71. 
T2. 


13. 


TH. 


15. 


76. 


17. 


78. 


Science of Explosives, edited by C. E. H. Bawn and G. Rotter, Part I, 


Ministry of Supply, London, 1956, (AD 127695). 


Torrey, Explosives Engineer, 1945, 23, Pp. 160. 


Huttl, Eng. Min, Jour., 1947, 17, p. 58. 


Clark, Amer. Inst. Min. Met. Eng., Tech. Pubs. 2157, 2158, 1947. 


Lawrence, Explosives Engineer, 1947, 25, Pe 171. 


C. F. Austin, "Lined Cavity Shaped Charges and Their Use in Rock and 
Earth Materials," Bulletin 59, New Mexico Inst. Min, and Tech., State 
Bureau of Mines and Mineral Research, 1959. 


C. F. Austin and J. K. Pringle, "Detailed Response of Some Rock Targets 
to Jets from Lined Cavity Charges," J. Petroleum Tech., January 1964, 


Cc. F. Austin, "Lined Cavity Shaped Charge and Its Use as a Drilling 
Tool," Trans. AIME, Vol. 220, p. 123, 1961. 


Draper, Hill and Agnew, U.S. Bur. of Min., R.I., 4371, 1948. 


Davidson and Westwater, Mine and Quarry Eng., May 1949. 


Michael J. Murphy, "Shaped Charge Penetration in Concrete: a Unified 
Approach," Doctor of Engineering Dissertation, University of California ~ 
Davis, January 1983, also Lawrence Livermore National Laboratory Report 
UCRL-53393, 1983. 


F. A. Bawn, R. P. Stanykovich, and B. I. Skekter, Physics of an 
Explosion, (AD 400151), New York, Research Information Service, 546, 1949. 


Charles E. Joachim, "Application of Shaped Charges to Rapid Atomic 
Demolition Munition (ADM) Emplacement," U.S. Army Engineer Waterways 
Experimental Station Miscellaneous Paper SL-83-20, December 1983. 


"Explosive Welding" Mechanical Engineering, May 1978. 


John S. Rinehart and John Pearson, Explosive Workin of Metals, Macmillan 
Company, NY, 1963. 


T. Z. Blazynski, (editor), Explosive Welding, Forming and Compaction, 
NY, 1983, 


Applied Science Publishers, 


VeoI5 


VI. The Generalized Equations 


The general conservation equations of mass, momentum and energy may be 
expressed as 


3 
yaa Ry, (eu); 


—J.2t. 
° De aX; (o> 


D 

a 
a 

a 
we 


the stress tensor equals the stress deviator tensor plus the hydrostatic 
stress. Also, 


E, = 1/2 Cu, u,) + EL or 


the total energy equals the kinetic energy plus the internal energy. 


The substantial derivative, or ,»represents 


Dele yy 

Dt at i ax, 
Thermal energy terms do not appear in the energy equation since for the 
micro-second time scale involved, there is no time to radiate or conduct heat. 


These equations constitute the basis of the hydrecodes. We need to 
relate the stress tensor to hydrostatic pressures and/or strains and strain 
rates. Then the strains and/or strain rates must be related to the velocities 
and/or the velocity gradients. The equations required for this step are known 
as the constitutive relationships which are intended to convert the stress 
tensor to a pressure and velocities and material properties (e.g., viscosity, 
shear modulus, etc.). These material properties are hopefully known, or 
expressible in terms of known quantities. 


Next, it is necessary to specify the appropriate initial and boundary 
conditions for the partial differential equations. It is also necessary to 
establish explosive burn routines and provide algorithms to couple the burn 
routine to the liner, for example, or to couple two complete sets of 
conservation equations such as for a penetrator and a target. Finally, a 
numerical solution technique and a stability criterion are required. This is 
obviously a complex procedure, but satisfactory results have been obtained 
from several available hydrocodes. Note that accurate metal and explosive 
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equations of state and accurate constitutive equations and associated material 
properties are prerequisite to a good solution. Assuming this information is 
available, a few limiting assumptions will be made. First, the following 
equations are available for, say, the behavior of a metal undergoing explosive 
loading. A method to treat the explosive is assumed. 


Equation Description Number of Equations Unknowns 
Continuity 1 Yom 1p, uy iz i, 2, 3 
Momentum 1 for each spatial direction 5-2, uss P i211, 2, 3 
Energy 1 5 =p, uy Ey isi, 2, 3 
Equation of State 1 3- 9 P, EL 


Thus, for a general, three dimensional problem there are six equations: 

1 continuity, 3 momentum, 1 energy and 1 equation of state. There are also 
six unknowns: 3 velocity components, the density, pressure and internal 
energy. Various simplifying assumptions change the situation as shown below: 


Case Equations Unknowns 

3-D 6 6 (op, P, Ep u,) iz, 2, 3 

2-D : 5 5 (p, Py Ey: u,) ies, 2 
(axisymmetric) (one less momentum equation) 
incompressible, 3 3 °(P, u,) iezt,2 
2-D (continuity and 2 momentum 

equations; p= constant) 

incompressible, 3 3 (CP, uy) isi, 2 
2-D, (continuity and 2 momentun No initial conditions 
Steady-state equations; = constant) required, one less 


independent variable. 


The axisymmetric case is the most commonly used, since this is usually a 
valid assumption for explosive-metal interactions. Note that a tremendous 
simplification is possible for an incompressible flow assumption. However, an 
incompressible flow implies a low Mach number flow and the complete absence of 
shock physics. This, of course, is not realistic. It may, though, provide 
approximate answers as will be seen. 
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Ye 


Equation Description Number of Equations Unknowns 


Continuity 1 Hoso, uy i211, 2, 3 

Momentum 1 for each spatial Sto, Uys P i214, 2, 3 Add P 
coordinate 

Energy 1 5 3 o> Uys Ey is 1, 2,3 Ada Er 

Equation of State 1 3: ps P, Ey Close System 


The point to be made here is that the generalized explosive-metal 
interaction is extremely complex. The governing equations are complicated; 
sufficient boundary and initial conditions must be specified; the constitutive 
relationships and material properties under the intense dynamic conditions 
encountered are not well known; equations of state must be accurately 
specified; and numerical algorithms and solution techniques must be accurately 
developed. 


Next, further simplications will be presented which lead to analytical 
models for the collapse of the shaped-charge liner. The next section starts 
with the one-dimensional hydrodynamic equation, i.e., the Bernoulli equation. 


Case Equations Unknowns 
3-D 6 6 (yo, P, EL» u,) 
i=1, 2, 3 

2-D 5 5 (oy P, Ey» u,) Need 1 less momentum 
(axisymmetric) i=, 2 equation. 
Incompressible, 3 3 (P, u,) p = constant. No energy 

2-D Continuity & iz, 2 or EOS needed. 

2 momentum. 

Incompressible, 3 3 No initial conditions 
2-D, steady-state As Above As Above required, one less 


independent variable. 


spatial coordinate 
time 
density 
velocity components 
stress tensor 
stress deviator 
Kronecker delta 
pressure 
total energy 


internal energy 


DEFINITION OF TERMS 
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VII. The Bernoulli Concept 


One of the most basic equations in shaped-charge jet collapse (and in jet 
penetration) is Bernoulli's equation. Due to the repeated application of this 
equation, and in order to clearly understand the assumptions inherent in 
Bernoulli's (and modified Bernoulli's) equation, it will be derived using 
several methods. The intent here is not to prove "Theorems" or even to 
introduce rigor in our methods, but to establish assumptions inherent to 
Bernoullits and similar equations. 


1. The Energy Equation 


Define a portion of fluid between sections (1) and (2) as shown below, 


@) 


u being the velocity. We assume one dimensional (1-D) flow and follow 
standard fluid texts such as Liepmann and Roshko [1]. 

This fluid portion can be considered to be bounded by pistons at 
(1) and (2) instead of a fluid. These pistons are equivalent to the fluid 
they replace, and help to clarify the work done on the system which consists 
of a definite mass of fluid. 


VIL=1 


During the small time interval in which the fluid is displaced to the 
position (1)! - (2)', a quantity of heat q may be added, (as external heat). 


Tne energy law yields, 
q + work done = increase in energy, 


and interpret the work done and the change (increase) in energy as; 
Work done: The volume displaced at (1) is the specific volume (per a unit 
mass) v,- For steady flow the displacement at (2) is also a unit mass, or 


specific volume, Vo. The work done by the pistons during this displacement is 


P, vy - P, ve, 


neglecting any external work between (1) and (2). 
Change in energy: The flowing fluid has internal energy e and kinetic 
2 
energy 1/2 u., per unit mass. After the displacement, there has been an 
increase of energy, e, + V/2 re corresponding to the displacement between (2) 
and (2)' and a decrease of energy, e, + V2 a,” from (1) to (1)'. 
Therefore, increase of energy = (e, + 2 us”) - (e, + 1/2 u,?) F 


For steady flow, the energy equation becomes 


2 2 

q+ Py Y- P, Va 2 en - e+ 1/2 (uy -4u, ) 

or for q = 0, the adiabatic energy equation yields 
P V2 ee P V2 é 

7%, +e + uy, =Pa¥ + ey + Us : 
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This equation relates conditions at two equilibrium states (1) and (2) of the 
flow. They are valid even if viscous stresses, heat transfer, or any 
nonequilibrium conditions exist between (1) and (2). However, (1) and (2) 
themselves must be equilibrium states. 


1f equilibrium exists all if equilibrium exists all along the flow (at every station), the 


equilibrium condition is valid continuously ar condition is valid continuously and 


Pv+er+ 1/2 a? = Const. 


Further, if the internal energy does not vary, 


Then 
2 
Pv + 1/2 u- = const., 
and if the density is constant (v = » ' 
P+ 1/2 6 u* = const., 


which is Bernoulli's equation derived from thermodynamic and conservation of 
energy considerations. ‘ 


2. Euler's Equation 
Now apply Newton's law to a flowing fluid, where 


Fema. 


Now view the flow in an Eulerian sense, i.e., observe the flow of a fluid 
particle as it encounters the various conditions in the tube through which it 
accelerates. 


The acceleration is in two parts: 


Since conditions vary along the tube, there is a velocity gradient in the 


direction of the flow, =. The rate of change of velocity is proportional to 


this gradient and the speed at which a particle moves through it. Thus, the 
acceleration due to convection through the velocity gradient is 


du 
ax 
Also, conditions at a given station may be changing due to nonsteady 
(nonstationary) effects or 


bu 
Che 
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Then, 
du du 
= fe Tyee 


a 
x ot ax 


Next, it is necessary to determine the force on a fluid element as shown below: 


0 
PA + 5, (PA) Ax 


PA — 


dab 


The force -2 (Ax A), and dividing by the volume of the element, A\x, and 


dividing through by the density p gives the force per unit mass or 
1 oP 


Roses 


x A Ox 


This result is valid for a particle (element) or any shape, 


¢ 
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converging, diverging or whatever. The proof is standard (e.g., [1]) using 
Gauss's law. Note that the force term neglects all viscous terms. Then, 


au gu _ 1 dp 
at * U Ox “5 dx 


which is Euler's equation. 


For steady flow, u we, dp = 0, since u =u (x), P = P (x), and for 
1-D, steady flow. ES 


For incompressible flow, (p = p_) and 


w + = const 
253 p = const., 


which is the classic derivation of Bernoulli's equation. Now consider the 
momentum equation and the flow through a control surface, where 


oO @ 


puas— (pualds 


puaA— — 
A 
peeigeead 
a 
pA a + pudA 2 = - A = (Euler's equation x pA) 


a a i as Pes ; : 
USE (eA) + u Ox (puA} = 0. (continuity equation x u) 


Adding these equations and with algebraic manipulation, we have 


a Din he Die aP BA 
at (puA) + a (ou"A) = - A rr oad 5x (PA) +P Poe 


for 1-D flow. 


Integrating between any two stations (1) and (2), yields 


2 
3 Z 2 as 2 
oC ; (puA) dx = oY, Ay ~ yy Ay = PAY - POA, + ; P dA , and 
2 
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The integral momentum equation is very general in that it is valid even when 
there are frictional forces and dissipation regions inside the control volume, 
provided they are absent at the reference stations 1 and 2. The integration 
of the differential momentum equation corresponds to a summation of the forces 
on adjacent fluid elements, and of the flow into them: 


Il! 
+ 


(pu’Alg 4 
2 {pu* A) (pu? A) 
The forces on adjacent internal faces are equal and opposite, and thus cancel. 
Similarly, the inflows and outflows through adjacent faces cancel, leaving 
only the force and flux at the control space boundaries. If there is a 
non-equilibrium region inside the control space, it does not effect the 
integral result. 


For steady flow in a constant area duct (A, = Ags 


which, sans the factor of 1/2, is similar to Bernoulli's equation. 


Harlow and Pracht [2] commented on the momentum conservation and the 
Bernoulli principle when considering jet penetration. They observed that the 
conservation of momentum principle, assuming the jet and target to be of the 
same width, was nearly the same equation as derived from the Bernoulli 
principle where the target is effectively infinite in width. 


3. Euler's Equation Extended 


Euler's equation can be found in any standard text, e.8-, Shapiro [3]. 
Euler's Equations of motion in 3-D cartesian coordinates (as derived above for 
a 1-D case) are in general: 


1 oP _ du ju au du _ Du 
Fay oe age a + Wot ar - (1) 
1 o3oP _ av av av av _ Dv 
poy ae ae YY oe eC (2) 
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and 


_ i sP _ aw aw ow aw _ DW 
Sora ae Get hae oe 


Consider 2 cases: 


1. Integration along a streamline (steady flow) [3]. A streamline is 
defined by the condition that at each instant the velocity vector is tangent 


to the streamline. Thus, 


‘ V 


Streamline 


dy sv dz wy 
dx u’dy y 
dx _u 

es aa] 


Now multiply equation (1) by dx, equation (2) by dy and equation (3) by dz, 


assume steady flow cm = 0) and for a particular streamline (using equations 


(4), 


du 


es 28 ax ey Bax + w Y ax iy 
ay Zz 


oe a OF eas 


from equation (1), and similiar expressions result from equations (2) and (3). 
Then, adding the resulting equations and noting that 


- Pays Pay 4 oP 
Pe 3g 8% + ay Wy 4 ay at 
and 

= ou 2Y gy 4 BY ae 
OU Sag Bh gy ee dz, etc 
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yields, 


- 2 dP = udu + vdv + wdw=d ( 


V being the resultant velocity. 
Thus, 
dP = -V¥ dV and for incompressible flow, . hs 


Vv 
P+P>” = constant along a streamline = Bernoulli constant of the 
streamline. ¥ 


This equation is only valid from point to point on a given streamline. Thus, 
the Bernoulli constant may vary from one streamline to another. Of course, if 
all the streamlines originate in a uniform flow region, the same Bernoulli 
constant would prevail over the entire flow field. 


2. Integration for Steady, Irrotational flow [3]. We abandon the 
restriction that only points along the same streamline are considered and 
assume only irrotational, steady flow. This means, 


av du 

a ay 8 8 

ow ov 

By” bz = 0 , and 
au aw _ 

ae ie 


Now, multiple equation (1) by dx, equation (2) by dy, and equation (3) by dz 
and use the irrotationality assumption,| which yields 


_ 1 oP _., ou du du _, ou av oY ay 
> ox dx =u Rx ox + Vv By dx + w BE dx =u Ox dx + Vv x dx + wsy ax, 


etc., and adding the equations, 
2 
eee eee tv oa a? 
ox 2 dea SD. 
etc., and finally 


2 
Th GES 
7 5 oP =d Ge 3 


which is identical to the Bernoulli equation obtained before. Thus, for 
steady, continuous, frictionless and irrotational flow the Bernoulli constant 
is the same for all streamlines and relates any two points in the flow field. 


The constant has a particular value for each streamline except for 
irrotational flow in which ease the constant is the same for all streamlines. 
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4, Navier-Stokes Equation 


We begin by investigating the derivation of the Navier-Stokes equations 
following, e.g., Schlichting [4]. 


These equations represent the equations of motion of a compressible, 
viscous, Newtonian fluid in three dimensional flow. The flow field will be 


specified by the velocity vector V = ui + vj + wk and the equation of motion 
(Newton's Second law) is 


Vv 
where pp is the substantial or Eulerian derivative of the velocity vector Vv. 


F represents the body forces (forces acting through the mass of the body) 
e.g-., gravity, and 


P represents the surface forces, stresses, shears, friction, and pressure 
forces. 


= 


We ignore F, 

The surface forces, P, depend on the strain rate of the fluid. The 
stress is related to the rate of strain via an isotropic, Newtonian fluid 
assumption. This means that the relation between the components of stress and 
the rate of strain is the same in all directions (isotropic) and this 
relationship is linear (Newtonian). However, shear stresses (viscous effects) 
will be ignored. 


To obtain the surface forces consider the fluid volume shown below. 


Vol. = dx dy dz 
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The x (plane 


one dimensional (1-D) case. 


_ oP 
5 -—X* where 
ax 
Poso i + * 
XOX Txy J 7 xz 


shear stresses (t). 


For a 1-D case and ignoring the shear, 


- A Oox 
Pel mq 
and neglecting body forces 
° Du _ 30, 
Dt ~ ox 


hydrostatic stress systems (fluid 
Stokes hypothesis will not be invoked, 
deformation will be neglected. 


to x) net surface force component is 3;_ 


oP 
‘dx dy dz for a 


The resultant surface force per unit volume is 


A 
k, after resolving P into normal stresses (.) and 


pressure) will not be considered, 
but deformations and rate of 


Thus, the normal stress will not be related to 


pressure, velocities or velocity gradients as done in the classic 


Navier-Stokes derivation. 


Instead consider a uniform stress 


field acting on the fluid element at 


rest (assume a negative a 0 somewhat analogous to a thermodynamic pressure. 


Then, expanding the substantial derivative and for incompressible flow 


o + S- ines 
at ax ax 
Du _ 3u bul 
for a l-D flow, (recall bot + ox 


Then, for a 1-D, steady state system 


o. 
ad ,u Xe. 
0 ak I tar 7 o 
or 
2 
u 
Pat oF const., 
with interpreted as a uniform normal 


ey ety y 
Ox! v ay Cr 
stress. This result is analogous to a 


"Bernoulli equation with strength" which with the basic Bernoulli equation is 
frequently cited in penetration and shaped-charge jet collapse and formation 


discussions. 


Of course, the full Navier-Stokes equations relax to Euler's 
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equations under the specified assumptions, which in turn lead to the classic 
Bernoulli equation as shown earlier. 


We conclude by remarking that the Bernoulli equation can represent a 
conservation of energy or a conservation of momentum depending on how it is 
derived. We also note that the "Bernoulli with strength" type equation, or 
the "modified Bernoulli equation" has some theoretical justification. 
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VIII. The Gurney Velocity Approximation 


The motion of a metal driven by an explosive was studied in the 1940's to 
early 1950's in a series of papers by Gurney, Sterne and Thomas [1,2,3,4,5] as 
applied to the motion of fragments resulting from the detonation of an 
adjacent high explosive [6]. This simple, approximate analysis assumes that 
the potential (or chemical) energy of the explosive charge before detonation 
is converted directly to the kinetic energy of the metal after detonation and 
to the expansion of the explosion products. The gas detonation products are 
assumed to expand uniformly and with constant density. 

The so-called Gurney approximation is based on the conservation of 
momentum and energy. The results represent excellent engineering 
approximations, within 10%, of the experimental results or detailed numerical 
results. The Gurney method may be applied to any one-dimensional 
explosive-metal interaction system. The amazing accuracy of the relatively 
simple Gurney formulae, over a wide range of metal mass (M) to explosive mass 
(C) ratios (0.1 < M/C < 10.0), is apparently due to offsetting errors. These 
errors are caused by ignoring rarefaction waves passing through the detonation 
product gases, which causes the calculated velocity to be too high, and 
assuming an initial constant density distribution of the detonation product 
gases rather than a distribution with a peak at the surface of the charge 
caused by the detonation wave, which causes the calculated velocity to be too 
low [6]. Based on references [6,7,8] the Gurney formulae will be derived. 


1. Open-Faced Sandwich 


From Kennedy (8*], a specific energy (energy per unit mass) 
characteristic of a given explosive is assumed to be converted from chemical 
energy in the initial state to kinetic energy in the final state. The final 
kinetic energy is partitioned between the driven metal and the detonation 
product gases by an assumed linear velocity profile in the gases. 


Figure 1 illustrates the open faced sandwich. It consists of a slab of 
explosive confined by a metal on one side only. This configuration is often 
used in experiments to obtain constitutive properties of materials or when it 
is necessary to impact a large target surface. The metal plate is termed the 
slapper or driver. The motion of the metal is obtained from momentum and 
energy conservation. For the open faced sandwich, Figure 2 illustrates the 
assumed liner velocity distribution. From this linear distribution, 

Y 
v AS (Y) = (Vv, + V) YS - Vv. Q) 

The y coordinate is interpreted to be Lagrangian or identified with a 
material particle. A given material particle is assumed to move at constant 
speed at all times for the purpose of constructing momentum and energy 
balances. The initial thickness and density of the explosive are Yo and Pe? 
respectively. The energy and momentum balances follow: 


Y 
2 
~1,,2 01 ° Y (2) 
CE = 5 Mv Fo, f [ie a1] dy , 


‘Note that this is J. E. Kennedy, not the Don Kennedy referenced earlier. 
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DETONATION 

PRODUCT GASES, 

METAL, M=TOTAL MASS 
'- PER UNIT AREA 


Figure 1. Open-Faced Sandwich. 
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VGAs, MAX = Y, 


EXPLOSIVE 
C= MASS/UNIT AREA 


Figure 2. Linear Velocity Distribution for the Open-Faced 
Sandwich Configuration [8]. 


VITI-3 


where 


Cs Toe 


E = Specific explosive energy 


mn 
t = metal thickness, and 
Y 
° y ; 
an eaveey f [« wy teva . 
° Y ? 
0 to) (3) 


and finally 
r 3 -1/2 
(1 +2 rae +1 M 
+a (4) 


6a + C , 


and 
os2tei (5) 


7 


Henry [6] provides additional mathematical details. Also, Henry presents a 
complete review of the Gurney formulae. The equations for some common 
symmetric and asymmetric configurations are presented in Figures 3, 4 and 5. 


Figure 3 depicts two asymmetric configurations, namely a flat sandwich and 
a cylinder. M denotes the metal mass and C denotes the explosive mass. These 
terms are obtained as the product of the density and thickness of the material 
in question. Thus, M and C are really masses _per unit area. The term V is 
the resultant metal velocity. The quantity Y2E occurs in all Gurney 
equations. It has the units of velocity and is termed the Gurney 
characteristic velocity for a given explosive. Gurney characteristic 
velocities have been tabulated for certain explosives by Kennedy [8], Henry 
(6], Jones et al [7], and many others, notably Dobratz [9]. Several 
explosives are tabulated from Reference [9] with their density, detonation 
velocity and Gurney characteristic velocity or Gurney constant. Note that two 
Gurney constants are given, one for warheads in which confining cases rupture 
at small expansions (prompt) and the other for warheads in which more ductile 
ease materials are used which expand further before rupturing (terminal). 


The tabulated densities are the actual density of the explosive used to 
determine the Gurney energy. The detonation velocities sometimes correspond 
to a different explosive density than that given in the table. It is 
recommended that the explosive of interest be accurately identified and the 
appropriate parameters be taken directly from Dobratz and the references 
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therein, or any other suitable source such as Meyer [10] or Reference [11]. 
La Rocea [12] presents a simple method for accurately predicting the Gurney 
constant and the relative power of any explosive. 


Figure 4 presents the spherically symmetric case and Figure 5 shows two 
asymmetric configurations, an open faced sandwich and an asymmetric sandwich 
which has two metal plates, one being a tamper. The derivation of these 
equations is relatively straightforward using the assumptions stated along 
with elementary calculus and algebra. Derivations are given in Kennedy [8], 
Jones, et al [7] and Henry [6]. Henry [6] considers a few additional 
explosive/metal geometries. The ratio of v//2E or the ratio of the metal 
velocity to the Gurney characteristic velocity is an explicit function of M/C 
or the metal to charge mass ratio. For the asymmetric sandwich where the 
second metal plate may be considered to be a tamper, v//2E is also a function 
of N/C, where N/C is the mass to charge ratio of the tamper plate. The 
velocity expression given is that of the metal plate, M. The velocity ratio 
v//2E is plotted as a function of M/C in Figure 6 for all the non-tampered 
eases presented. Since values of VY2E have been tabulated for certain 
explosives, estimates of the metal velocity can be obtained for a given M/C 
for the geometries considered here. 


Figure 7 plots the metal velocity as a function of M/C for an open faced 
sandwich and a few values of the Gurney constant, Vee, corresponding to 
PBX-9404, TNT and Comp B. This plot allows easy interpolation or 
extrapolation to other values of the Gurney constant. 


The range of applicability of the Gurney formulae is restricted due to 
the simplifying assumptions in the derivation. These restrictions, according 
to Kennedy [8], Jones, et al [7] and Henry [6] are listed and discussed in the 
following tables. The Gurney assumption of a linear velocity profile of 
constant density detonation product gases greatly deviates from eonventional 
gas dynamic theory. This assumption introduces the largest error in 
configurations involving a free explosive surface such as the open faced 
sandwich, for which the Gurney method may overestimate the metal velocity and 
impulse. Henry [6] used a parabolic pressure profile in the detonation 
products to remove the assumption that the gas density was constant at any 
given time. A simple gas equation of state was also used. Henry concluded 
that the additional mathematical complexity, inconsistent with the basic 
simplicity of the Gurney method, was not worthwhile based on the small 
improvement obtained in accuracy. Other extensions of the Gurney method are 
given in Reference [13]. 
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Explosive Density D veE, km/s 


gn/om> kn/s , Prompt Terminal 
5-7 om 19 - 26 mm 
Comp A-3 1.61 8.47 2,402 
1.59 2.63 
Com B 1.71 2.70 
1.717 7.89 2.35 2.756 - 2.821 
1.717 2.71 
Cast 1.68 2.402 
4.62 2.32 
Pressed 1.59 2.335 
Comp C-4 1.52 8.37 2.176 
Comp C-3 1.60 7.63 2.68 
Cyclotol 75/25 1.754 8.2 - 8.3 
Cast 1.69 2.286 
Pressed 1.64 2.362 
DATB 1.68 7.52 1.975 
Explosive D 1.50 6.85 1.942 
HBX~1 1.70 7.31 2.213 
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Explosive 


HBX=3 


RMX 


LX=-14 


NM 


NQ 


Octol 75/25 


1.81 


1.89 


1.835 


1.14 


1.44 


1.81 


1.821 


DETASHEET C (DuPont)1.48 


(BRL) 


km/s 
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"BE, km/s 


Prompt Terminal 
5 - 7 mm 19 = 26 mm 
1.984 
2.97 
2.80 
2.41 
1.896 
2.58 
2.83 
2.1 - 2.3 
2.8 


Explosive 


PBX-9011 


PBX=9404 


PBX-9502 


Pentolite 50/50 
Cast 


Pressed 


PETN 


RDX 


TACOT 


Tetryl 


TNT 


Cast 


Pressed 


Density 


8.5 


8.8 


7.71 


6.73 
6.93 
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Y2E 


Prompt 


5-7 om 


2.82 


2.301 
2.317 


2.451 


22274 


2.039 


2.097 


2.103 


, kn/s 


Terminal 


19 - 26 mm 


2.50 


2.419 - 2.505 
2.37 
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Restriction 


Range of M/C ratio 


Acceleration Phase 


Direction of Detonation 
Propagation 


Gas Velocity Profile 
Assumption 


RESTRICTIONS ON THE GURNEY MODEL 


Remarks and Recommendations 


Henry [6] claims accuracy for 0.1 < M/C < 5 and 
Kennedy [8] states a range of 0.2 < M/C < 10 for 
velocity calculations. Impulse calculations are 
usually acceptable for M/C > 0.2. 


The Gurney method in its basic form is not capable 
of analyzing motion during acceleration. However, 
Henry [6] and Jones, et al [7], used the 
assumption of a liner velocity profile in the 
uniformly dense gases in conjunction with an 
equation of state of the gases to yield 
acceleration solutions for flying metal plates. 
The detonation product gases must be allowed to 
expand sufficiently to complete the acceleration. 
The flying metal will reach its calculated 
velocity only if no external forces (or 
interactions) are applied during the initial 
acceleration phase. 


The detonation of the HE drives the metal at a 
given velocity (approximately) for a given M/C 
regardless of the angle between the detonation 
front and the metal. The direction in which the 
metal is driven will vary slightly with the angle 
(see Taylor angle discussion). 


The assumed linear velocity profile and constant 
density gas expansion are major assumptions. 
Ignoring effects of both rarefaction waves and of 
pressure peaks near the metal surface would seem 
to be cancelling errors. However, these 
assumptions are inherent to the simplicity of the 
Gurney method. 
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RESTRICTIONS ON THE GURNEY MODEL 


Restriction Remarks and Recommendations 
One-Dimensional Motion The Gurney approach is not valid for the estimate 


of variations in local velocity of a plate driven 
by a charge with a tapered thickness (or a tapered 
metal thickness). In this case, an average value 
of M/C could be used to estimate the final 
velocity of the entire plate. 


Metal Strength Effects Forces exerted by the metal to oppose deformation 
are not considered, other than inertia. Hoop 
stresses in cylinders and spheres can reduce the 
metal velocity in explosions and the strength 
effect is greater in implosions. 


Metal Spallations Metal spallation may occur when M/C < 2 for high 
density explosives and metals. Spallation ean, in 
come cases, be avoided by introducing an air gap 
of a few millimeters between the explosive and the 
metal. This results in a small decrease in metal 
velocity. 


Early Case Fracture The leakage of the detonation product gases 


through fractures in the metal case can decrease 
the final metal velocity perhaps significantly 


depending on the nature and extent of the leakage. 
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The specific impulse delivered by the explosive can be calculated as the 
total momentum imparted to the metal body divided by the total explosive 
charge mass or 

ag 
year “7 (6) 


The specific impulse, Typ) will now be derived for an unconfined surface 


charge. Explosive may be detonated directly on the surface of.a massive metal 
body in order to deliver a desired impulse for testing pruposes. If the 
loaded body is considered to be rigid, all detonation product gases will flow 
away from the surface and a maximum specific impulse will be delivered to the 
body. For a very large metal body mass to explosive mass ratio, M/C >> 1, the 
open faced sandwich velocity, equation (4), reduces to 


3.6 
Vu V2E IM ’ (7) 
3E 
then tae V2 , from equation (6). 

Also since 
Vy =v [: - + | from equation (5) and for M/C >> 1, 

Vo = 2V M/C and with V given by equation (7), 

Vio = 2 V3 V2E or 

° 4 

v= V6E 

e (8) 


for the maximum free gas velocity. 


The impulse delivered to a body by the explosive can be increased by 
tamping or confining the explosive with an outer layer of metal. Tamping acts 
to hold back the detonation gases, keep the pressures high and force the 
detonation gases to do useful work against the body. A tamping and a 
confinement are analogous. The Gurney method can be used to estimate the 
change in impulse as a function of the tamper areal density. If the heavy 
metal body sandwich is considered to be rigid, it resembles the plane of 
symmetry in the symmetrical flat sandwich configuration shown in Figure 3. 

The velocity of the tamper plate is then given by the flat sandwich of Figure 
3 with the tamper to explosive mass ratio N/C, used instead of M/C or 


- 1/2 
vee [Eo 3] t . (9) 
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The effective specific impulse of the explosive for a tamping ratio of N/C 
follows from the equation for the asymmetric sandwich of Figure 5, namely, 


(10) 


Assuming M/C >> 1 and with A = : - 5 we ; the specific impulse becomes 
1/2 
16? tamped = /2E (N/C + 1/2) / (N/C + 1/3) ; (11) 


Note that when N/C = 0, i.e., no tamper, equation (11) reduces to 


3 
I w= as derived earlier. 
sp 


The acceleration of explosively driven metal plates can also be 
determined by using the Gurney assumption of a liner velocity profile in 
uniformly dense gases in conjunction with an equation of state for the gases. 
Acceleration calculations showing the displacement as a function of time or 
speed are important when estimating the effects of venting in a perforated 
casing. Henry [6] and Jones, et al (7], calculate the acceleration of the 
casing for various sandwich configurations. These same authors also extended 
the Gurney analysis by considering the pressure gradients within the gas in 
order to remove the assumption that the gas density is constant at any given 
time. These results however, complicate the simple Gurney model. 


The Gurney model yields explicit algebraic relationships for estimating 
the velocity imparted to a metal in contact with a detonating explosive. The 
model can be used for several simple geometric metal/explosive configurations. 
The specific impulse of an explosive oan also be obtained for these geometries 
and is directly related to the Gurney energy of the explosive. Thus, the 
Gurney method can be directly and simply applied to a multitude of 
explosive/metal interaction problems. Design and parametric studies can be 
readily performed. Applications include warhead and fragmentation design, the 
study of the efficiency of conversion of the high explosive chemical energy to 
the kinetic energy of the plate, explosive initiation by the impact of an 
explosively driven plate, the launching of a dielectric plate by electrical 
detonation of a metal foil, and the calculation of layers of metal fragments 
in conjunction with shock wave physics. 


Note that the Gurney method does not consider the propagation of shock 
waves through metals, but only the terminal effect of the shock propagation. 
Duvall [14] provides an excellent description of shock wave propagation, 
transmission and reflection. He describes a plane wave lens, a plane wave 
generator and methods of forming oblique shock waves. He also states the 
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pressure levels attained by the various methods of shock wave generation, 
Duvall's article [14] explains shock wave phenomena without recourse to 
advanced mathematical concepts. Thus, Reference [14] is an excellent 
introductory article for the layman. Courant and Friedrichs [15] provide an 
excellent mathematical treatment of shock wave physics. Shock wave 
Propagation through solid media is an important aspect of explosive/metal 
interaction, but is beyond the scope of this introductory report. 


An extension of the Gurney method was given recently by Chanteret [16]. 
Chanteret developed an analytical model for symmetrical geometries with energy 
partitioning during expansion. This model allows a calculation of the Gurney 


energy as a function of the expansion ratio. The expansion ratio is (R/R,)" 


where R is the current value of the radius of the explosive-metal interface 
and Ry is the undetonated radius of the explosive-metal interface. The 


constant n depends on the geometric configuration (n = 1 for symmetrical plane 
sandwiches and n = 2 for cylinders with a constant wall thickness). 
Chanteret's results can be used in conjunction with the classical Gurney 
formulae to obtain the entire velocity-time curve for all symmetrical 
Steady-state configurations. Also Chanteret extended the Gurney model to 
include imploding cylinders by introducing a fictitous rigid boundary in the 
explosive cylinder. The results are claimed to agree well with available 
experimental data and two-dimensional hydrocode results. 


The Gurney method and extensions of this method are useful analytical 
tools. The flyer plate velocity depends only on the M/C ratio and the Gurney 
constant or Gurney characteristic velocity, /2E. This constant is available 
for several explosives (as we have seen), but when the Gurney constant is not 
available, Kennedy [8] recommends that E ~ 0.7 Hp be used when Hp is the heat 


of detonation. This recommendation follows from the fact that 0.61 < EH < 


0.76 for most explosives. Another option is to use La Rocea's method [12] to 
calculate the Gurney constant. 


The Gurney method ignores rarefaction effects and the shock 
characteristics of metals. Thus, we can expect errors to occur when these 
effects are important. An example would be a short open ended eylinder which 

‘will deviate from the Gurney prediction near its ends. Also, the difference 
in shock impedance between different metals will alter the energy coupling and 
as a result, the same M/C ratios with different metals but with the same 
explosive will yield different velocities. The applicable range of M/C values 
is 0.2 < M/C < 10 as discussed earlier. Note that as M/C decreases, e.g., 
when adding HE to a given metal slab, the velocity cannot increase 
indefinitely as predicted by the Gurney method. 


Finally, the Gurney equations assume that the metal moves in a direction 
normal to its surface. This is true when the detonation wave encounters the 
metal at normal incidence. It is not true when the detonation wave encounters 
the metal plate at a grazing incidence. The grazing incidence is treated by 
another model known as the Taylor angle approximation. 


Before investigating the Taylor angle approximation, a few Gurney type 
formulae relating to flat plates and imploding devices will be discussed. 


VIII-18 


For shaped charges, a popular formula for a single flat plate backed by 
slab of explosive is given by [17], where = M/C, and ve is the liner 
collapse velocity, and 


_ F 1/2 
V, = 22 | ————_. 
re ee (12) 


For this same configuration, with the detonation wave (of velocity D) 
propagating in a direction tangent to the liner, Duvall and Erkman [18, 19] 
used hydrodynamic theory to obtain the collapse velocity formula 


Vv, 22 [: : =~ ( iWies rea | : (13) 


Kleinhanss [20] also presented equation (13) and credited Trinks with another 
formula based on semi-empirical data related to the explosive launching of 
plates, 


ce 
ue 0.36D arctan (SF - (14) 


These formulae were developed for flat plates, and although they are 
quite accurate, they do not account for the effect of curvature such as in 
implosive geometries, e.g., a conical shaped charge. Kleinhanss [20], in 
experiments with imploding cylindrical charges, showed a strong dependence of 
the liner collapse velocity (Vv) on its radius. He obtained the 
semi-empirical formula 


Se.3 e(2z, - e) rT 
; ee 1 
fe) (Cy + € 2(d)) (15) 


where ¢ is the metal thickness, b = roo hye? and ry being the outer and 


° 
inner explosive radii, respectively, and thus b is the explosive thickness. 
cy and f(b) are empirical constants that depend on the type of explosive and 


metal being used. Other studies on implosive geometries were given by 
Chanteret [16] as discussed earlier, Chou et al, [21] and Hennequin [22]. All 
of these studies provide improved results over the basic Gurney model for 
implosive geometries. 


Other formula are available for plate acceleration and are often used in 
shaped charge liner collapse and explosive welding models. For example, Singh 
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(23] considers the velocity of the casing of a thin walled cylinder filled 
with explosive and detonated at one end. The velocity of the casing at the 
time of fracture is given by 


(16) 
where each element of the casing subtends an angle q with the axis of 
propagation of detonation. 


Deribas [24] gives 


where 1 is the charge to mass ratio. 


Shushko, et al [25], gives 


Vv 2 on ag 
ee se Ne eet aa 
n(n + 6) 


Mikhailov and Dremin [26] provide several frequently used expressions for 
the plate speed. They are: 


V iLo+ zen - 1 
“o LY 27 a (19) 
D a ; 
qa? 
0.602, ana 
o 2e+n ? (20) 
D oe 5 BES (21) 
161 a 5 ie 
where y = WT O= 14+ 2¥)/C1 - Dt » t is the time taken to 


travel the flight distance or the length of the driver plate, and H is the 
height of the explosive charge. 
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Smith, et al [27] used formulae analogous to equations (12) and (17) to 
determine the terminal velocity of flyer plates used in plate acceleration 
experiments. Smith, et al [27] also calculated the flyer plate terminal bend 
angle which will be discussed shortly. Both the bend angles and flyer plate 
terminal velocities are in agreement with the experimental values obtained for 
three different explosives. 


Explosive-metal interaction models are also useful in fragmentation 
calculations. Fragmentation, of course, does not always result from an 
explosive-metal interaction. For example, fragments are generated when gas 
pressurized pipes or vessels fail. The velocity of the fragments following a 
gas pressurized system failure is necessary in order to determine the damage 
to neighboring equipment on structures. Fragment velocities resulting from a 
pressurized pipe or vessel failure are usually predicted by the Moore equation 
[28]. The Moore equation is semi-empirical based on experimental data from 
the velocity of explosively generated fragments (i. e., Gurney data). 


Another method for estimating the fragmentation velocity from ruptured 
pipes or vessels follows from the theoretical model of Baum [29]. 


Of course, other Gurney like models are also available. The section on 
explosive welding and the sections on jet formation provide references to 
other metal acceleration models, Walters [30] and Walters and Harrison [31] 
provide several references (mostly Soviet), related to the acceleration of 
metals by high explosives. Dehn [32], recently published a report containing 
an overview of Gurney type velocity approximation models as well as the Taylor 
angle (grazing incidence) models. Dehn provides additional historical 
information and correctly analyzes geometric configurations in this report 
such as the Jelly Roll, the Dagwood and similar explosive-metal multilayer 
arrangements. 


In conclusion, care must be taken in comparing the various 
explosive-metal interaction formulae. Many authors use M/C (ratio of the mass 
of metal to the explosive mass) as a key parameter. Other authors prefer C/M 
as the parameter. 


Next, grazing detonation waves will be investigated, 
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THE TAYLOR ANGLE APPROXIMATION 


The Taylor approximation [33] 1s shown in Figure 8, which depicts a 
grazing (parallel) incidence of the detonation wave to the metal surface. The 
plate is deflected at an angle 6 from its initial position. Acceleration of 
the metal to its final velocity is assumed to be instantaneous. For this 
steady-state condition, the metal plate is assumed to undergo pure rotation, 
i. e., no net shear flow, or no change in length or thickness. Thus, the 
plate element that was at P initially will be a P' after launch and the 
lengths 


P = OP! 


Construct a line from 0 perpendicular to PP'. This line bisects the angle 
since OPP' is an isosceles triangle. If the time, t, is measured from the 
time when the detonation wave passes point P, then: 


OP = Dt, 
PP! = Vt 
and 0 _ PP'/2 _ Ve JV (1) 
siny = So en fs 
2 OP 2Dt 2D 


An estimate of the direction in which the metal plate is projected is the 
angle 9/2 which can be determined from equation (1) with V known from the 
Gurney method and D determined from the specified explosive. As the plate 
moves, it will be tilted an an angleé from its initial position, and by 
assumption, will not be rotating. 


When using a streak camera located perpendicular to the initial charge 
axis, a component of the velocity, namely Va, the apparent velocity, was 
measured by Kury, et al [34]. This velocity can be related to V by the 
geometry or 


Va = D tan®, (2) 
then 


V__ 2D sin 6/2 _ 2_sin 6/2 _ cos 6 (2 sin 6/2) 


Va D tan 6 sin 8 2 sin 5 cos cs 
cos 8 2 2 


(3) 


Vv cos 6 


or Va cos 6/2 


Hoskin, et al [13] expressed their results in terms of the velocity 
component Ve the velocity component normal to the flight plane of the plate, 
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since this velocity is of interest when the flying plate is used in impact 
studies. Since the angle between va and V, is @, we can write 


N 

Vn = D sin 6, and (4) 
V__ 2D sin 6/2 _ 1 7 

V, (bsince > * eoseye wee O02 i (5) 


The difference between V, Vy and vy is usually only a few percent. Also, the 


term V/2D is approximately the same for many explosives. Thus, 6 is 
approximately constant. If we fix the mass of the plate (its thickness) being 
driven, we observe that increasing the detonation velocity of the explosive 
implies an increase in the plate velocity such that V/2D is approximately 
constant. These concepts will be used to study jet formation theory. 
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IX. JET FORMATION 


Birkhoff, Mac Dougall, Pugh and Taylor [1] formulated the first theory of 
conical shaped charge jet formation assuming that the detonation wave produces 
such large pressures during the collapse process that the material strength of 
the liner may be neglected. In fact, the liner is treated as an inviscid, 
incompressible fluid. The conical liner is modeled as a wedge, and a steady 
state collapse model is assumed. Thus, the liner elements are instantaneously 
accelerated to their final collapse velocity. The steady state model 
predicted a jet with a constant length equal to the slant height of the cone. 
It has been observed however, that shaped charge jets possess a velocity 
gradient with the tip traveling much faster than the tail, causing the jet to 
stretch and eventually break up. This steady-state theory was later modified 
by Pugh, Eichelberger and Rostoker {2] to include the jet velocity gradient. 
The modified "non-steady" theory is based on the same principles as the 
original theory except that the velocities at which the various liner elements 
collapse is not constant, but depends on the original position of the element 
in the liner. We begin with an examination of the Birkhoff, et al [1] steady 
State model. 


As the detonation wave impacts the liner the pressure on all sides of 
the liner is assumed to be equal and the liner walls are assumed to collapse 
inward at a constant velocity, say Vo: The angle 28 between the moving walls 


is greater than the original apex angle 2c, because of the finite time 
required for the detonation wave to sweep the liner surface from apex to base. 
Figure 1 describes the geometry of the collapse process. The conical liner is 
symmetric so a represents one half of the liner apex angle and g represents 
the collapse angle. The liner assumes a velocity, Vy which bisects the angle 


APP' of Figure 1. To show this consider a coordinate system having a constant 
velocity such that the origin moves from point P to Pt in a unit time. In 
this coordinate system steady state conditions exist at the origin and the 
liner moves inward along P'P and flows outward along the path PA. The 
velocity of the liner passing through this region changes its direction but 
not its magnitude, because the pressure forces are everywhere perpendicular to 
the motion. P!'P and P'B represent the entering and emerging velocities, 
respectively, of the liner in the moving coordinate system. P'B is 
constructed parallel to PA and P'P and P'B are equal in magnitude. Since the 
velocity in the moving coordinate system is PP' the velocity of the collapsing 
liner in the stationary system is the vector sum , 


PP' + P'B = PBsV_. 


Since P'P = P'B, the triangle BPP' is isosceles and since P'B is parallel 
to PA, angle BPP' = angle PBP' = angle BPA. Call these angles 6. Thus, Vo 
bisects angle APP’. 

The walls of the collapsing liner are two planes moving inward. The 
junction of these planes moves from A to B with a velocity Vy 


B- & 
2 


) /sin B. (1) 
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From Figure 1, the angle OPA = ¢ follows from the fact that o + ¢@ + 180 - g = 


180 or ¢ = 8 - Oo, i.e., the sum of the angles of a triangle equals 180°, 
Next, consider the triangle PXB where 180 = 90 + 6 + ({XBP = y) or y = 90 = 6. 
In triangle XAB, 180 = 90 + 8 + XABX or JABX = 90 - 8. Then along the line 


OPP', 26 + ¢ = 180 or 6 = 90 = $/2 = 90 -(£ 5 &“ \, from above. Then since 
y+ 6= 90, y+ 90 ~-(£52)s 90 or y= = _ Angle PP'B = $ since in 


B - 
2 


triangle PP'B we have 26 + YPP'B = 180 or 180 - 6+ 4PP'B = 180 or YPP'B = 9, 
This establishes the geometry of Figure 1. Now for some trigonometry. 


In triangle APB by the law of sines, V,/sin a= V,/sin fp or vy = Vy sin 6/ 


sin 8 The sin 6 = sin (90 -(254 ) = cos (8 - a ) and 
2 


(2) 


V2 =V, cos (g 3 ©) ei B= V, cos (/2)/sin B. 


A moving observer positioned at A would observe any point P in the upper 
plane approaching him with a velocity v, cos g++ Vy cos g. This velocity, Vor 


becomes V, = V) cos B+ V, cos (50 - (2 = «)). V, cos g+ Vy sin (58), and 
from equation le), 4 


= 8 - a : 
Va, = v.. [oe ( 5 ) /tan 8 + sin (3) 
Also Up = U cos g where U is along PP' and using the law of sines for triangle 
v 
PBP', ° -—_¥U_ or 


Sin (6 - a) Sin 6 


U* Sin (@- a) 5" 20 - ( 7)" ~sin =a)” 
then 
U ux ( S22) 
Do. 28 2 
cos a Ssin(B - a) . (4) 
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Returning to V3 from equation (3) the observer sees a "jet" moving to his 


right and a "slug" moving to his left, as in Figure 2. Also, as viewed by 
this observer, the whole process appears to be unchanged by the lapse of time, 
or steady state motion exists. This, with the earlier assumptions of 
inviscid, incompressible and one-dimensional flow allows the use of the 
Bernoulli equation or 


P=%o ur = constant, (5) 


which relates the pressure and the corresponding velocity U. The pressure at 
any point in the liner determines the velocity at that point. Assume the 
liner moves away from the detonation so fast that the pressure of the surface 
‘decreases rapidly and the pressure on all surfaces of the collapsing liner is 
constant. This allows free streamlines, i.e., the boundary streamlines are at 
constant pressure and density (and hence velocity). As viewed by the 
observer, the jet and slug will appear to recede with exactly the same speed, 
Vos as the walls approach as in Figure 1. 


Returning to the stationary system of coordinates, it is seen that the 
jet, traveling to the right in Figure 2, has a velocity 


VeV,+/V,, 

whiid the“slug, traveling to the left in the moving system of Figure 2, 
has a velocity to the right given by 

ve = v, - vo: 

To further visualize this process consider that the point P (fixed in the 
upper plane) travels to point B (fixed in space) in unit time, the material 
from the inner surface of the upper plane included between PA and AB moves 
into the jet, and the front of the jet moves to the right a distance PA + AB 
in the same time, forming a high velocity jet. This jet velocity is 
V = V., + V., and from equations (2) and (3): 


1 2 
E - a) B= «| 
cos \—3— cos\—— ere (2 ong ) (6) 
VEN sin B * “tan 8 2 : 


Also, the outer surface of each plane forms a slug moving with the lower 
velocity, vy, = v, - V5 or 


x \ 
Q- B- a. 
cos (e] ets) ets (2 - a ) } (7) 
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Conservation of momentum dictates the division of material between the jet and 
the slug. Let m be the total liner mass per unit length of the two collapsing 
planes (sides of the wedge) approaching the junction. Let my represent that 


part of the liner mass entering the jet and nm, represent that part entering 
the slug, so 


mem, +m, . (8) 


Equating the horizontal momentum components entering and leaving the 
junction A of Figure 2 in the moving coordinate system yields 


V = _ 
mV, cos 6 nm, Vv; ms Vy : (9) 
Solving equations (8) and (9) simultaneously gives 


ue = (m/2) (1 - cos g) , 


(10) 
and 


mo = (m/2) (1 + cos B) . (11) 


According to this model, the velocities of the jet and slug and their 
cross-sectional areas are constant. 


The analysis assumed a wedge configuration. A conical liner may be 
treated in the same way. In the conical case, the walls converge on the axis 
from all sides. The moving observer must travel at the same rate as in the 
case of the wedge. In order for the process to appear stationary to him, the 
total mass per unit distance along the axis must be constant. This is 
approximately true for a constant wall thickness liner. In reality, the liner 
wall thickness would have to be tapered to be inversely proportional to the 
distance from the apex. 


For the detonation wave traveling parallel to the axis with speed Uns as 
in Figure 1, equations (6) and (7) become, from equation (4), 


2 (25+) : 
V= Zosa San (8 - a) | csc B + ctn B + tan 5 > (12) 
and 
U 
2 (2:2) 
2 - t ’ 
Mg ora a a 2 (13) 


The jet velocity, V, increases as q decreases, since 8 also decreases. The 
jet velocity approaches a maximum as 9 > 0, or V = Up {1+ cos + tan 8/2) or 


68> 0asarO0soV=z2 Up» which means that the jet velocity cannot exceed 
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twice the detonation velocity. Also, as 27 8-70, Vg 0. Note that as 


@ > Q the conical liner approaches a cylinder. Cylindrical lines are capable 
of generating high velocity (and low mass) jets. 


For the hypothetical case of a conical wave front moving perpendicular to 
the surface of a conical liner such that the wave strikes all surfaces at the 
same time, 8 = and the velocities of the jet and slug from equations (6) and 
(7) this time, become: 


Vv Vv 
0 oO 
= —2_ ¢ = 2 - 
Vv sina (1 + cos aj, and Vo sino (1 - cos a} 


With this type of wave front the jet velocity could be increased indefinitely 
by decreasing o. However, as& > 0, Vo > Q, my > 0 and the jet momentum 
mV 
n, Ve > 
Finally, the theory of Birkhoff, et al predicts steady state jet and slug 
velocities as in equations (6) and (7) or equations (12) and (13) and masses 
as in equations (10) and (11) for conical or wedge shaped liners. 


sina> 0, 


Comments on the Birkhoff, et al Solution 


The steady state model provides quantitative agreement with flash 
radiograph experiments. The model tends to overpredict the jet velocities. 
Also, the jet possesses a velocity gradient from tip to tail and continues to 
stretch after the walls have collapsed giving a jet length greater than the 
slant helght of the cone. The ductility of the jet under the intense 
pressures and dynamic collapse conditions is responsible, at least in part, 
for the elongation of the jet. 


Important modifications have been made to the steady state theory of 
Birkhoff, et al notably by Pugh, Eichelberger and Rostoker [2] and Godunov, 
Deribas and Mali [3]. Pugh, et al, hereafter referred to as the PER theory, 
developed a non-steady theory based on the same concepts as Birkhoff, et al 
[1] except that the collapse velocities of the various liner elements are not 
the same for all elements but vary depending on the original position of the 
liner element. Thus, the collapse velocity decreases continuously from the 
cone apex to its base, producing significant jet elongation. 


In the Soviet Union Godunov, et al [3] hereafter referred to as the 
visco-plastic theory, modified the steady state theory to include 
visco-plastic behavior, i.e., strain rate dependence. 


Both of these models will be discussed in turn. 


IX-7 


The PER Theory 


The PER theory assumes a variable instead of a constant collapse velocity 
for the walls of the conical (or wedge) liner. This condition greatly 
improves the steady state theory results. The collapse velocity decreases 
from the apex to the base of the cone and Figure 3 illustrates the effect of 
these velocity variations. Note that as the collapse angle 8 increases, the 
jet velocity decreases, but the portion of the liner entering the jet 
inereases. Figure 3 also shows how the decreasing collapse velocity increases 
8. As the detonation wave travels from P to Q along the conical surface APQ, 
the liner element originally at P collapses to J. The liner element 
originally at P', starts later and collapses slower than P, and arrives at M 
at the same time P arrives at J. The element P' would have reached N when P 
reached J if their collapse velocities were identical. Thus, with a constant 
collapse velocity, the surfaces remains conical, or QNJ is a straight line. 
However, since P' has a slower collapse velocity than P, the collapsing liner 
has the non-conical contour QMJ of Figure 3. The angle 8 is greater than Bt, 
the "steady-state" 8. This assumes that each liner element is thin and not 
affected by its neighbors, consistent with the hydrodynamic assumption. 


Next, Figure 4, which is similar to the steady state model of Birkhoff, 
et al, illustrates the flow situation. QJ is parallel to PA and equal in 
length to PQ. If the magnitudes of QP and QJ are equal to U, they represent 
the velocities in the moving coordinate system of the liner elements entering 


and leaving the region P. The vector PJ = Vy is the velocity of the 


collapsing liner element in a stationary system of coordinates. Thus, the 
arguments of Birkhoff, et al [1] are valid. The liner element does not move 
perpendicular to its original surface, but along a line that makes a small 
angle 6 with the normal. This angle, from Figure 4, is 


sin 6 = V,/2u (14) 


If Vo is constant, B = 8* and 6 = (£54) and all the PER results relax to 


the steady state model. 


By the proper choice of coordinate systems, the geometric relationships 
at the moving junction (J) are shown in Figure 5. The axis of the cone is 
along JR and OJ is the element of the liner moving toward the axis. This 


element has a velocity OJ Vin moving coordinates. The velocity of the 


moving coordinate is JR = V5. Using the law of sines from Figure 5 we have 
Ve Vy cos (a + 6)/sin 8, and (15) 


Vv. = V_ cos (B - « - 6)/sin B. 
1 ° (16) 
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Figure 3. The rye Process for a Variable Collapse Velocity Liner, 
PER [2]. 
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Figure 4, Velocity Vectors of a Collapsing Conical Liner Element, PER [2]. 
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Figure 5. Relationship Between Vos The Liner Collapse Velocity, v,The 
Collapse Velocity Relative To The Collision Point, and Vy The 
Collision Point Velocity. 
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In fixed coordinates, the velocities of the jet and slug are yy = y, + V and 


ve = v, - V, respectively. Using equations (15) and (16), along with some 


trigonometric manipulations, the jet and slug velocities become: 


4 (17) 
J 


Vv 
8 


Vy (csc 8/2) cos (a + 6 - 6/2) and 


v (sec 8/2) sin (a + 6 = 8/2) , (18) 


respectively. Note that for B = 6* sa + 26 we get the Birkhoff, et al jet 
and slug velocities, after some trigonometric manipulation. 


We can eliminate from equations (17) and (18) using equation (14) to 
yield 


(19) 
vi = V, (ese 8/2) cos (a - 8/2 + sin! Vj/2u) , 


and Ve = Vo (sec 8/2) sin (a - 8/2 + sin! V,/2u) 
(20) 


These equations are valid in either the steady-state case when Vy is constant 
or in the nonsteady case where Vo varies. In the steady state case, however, 
f could be expressed in terms of a, U and Vy and need not appear in the 
equation. 


Of the four unknowns Vis Vor m, and m, we still need expressions for m 


J 


and n.. These equations follow from the conservation of mass and momentum, 


or, where m is the liner mass and ms and my are the jet and slug masses 
respectively, 


dm = dm, + dm, 
8 


j 
dm. 
— sin’ (8/2), and (21) 
dm 
Ss _ 2 a (22) 
in = cos” (8/2) . = 


Equations (21) and (22) are identical to the steady state theory, (equations 
(10) and (11)}. 


Equations (19), (20), (21) and (22) describe the mass separation and the 
velocities of each element of a conical liner. They depend upon the cone 
angle 2q.and the detonation velocity Up = U cosa, and on the collapse angle 8 


and velocity Vo» both of which are variables. 
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From Figure 3, observe that the calculation of the angle Be would be more 
straight forward than the calculation of g, since Vo is different for each 


element of the liner and thus the contour of the collapsing liner is not 
straight but curved as line QMJ of Figure 3. Toward this end, let the 
cylindrical coordinates of M be (r, Z) and the coordinates of P' be (X tan %, 
X) in Figure 3. The coordinate directions are indicated in Figure 6. Then we 
have 


ZeX+V, (t - T) sin A, and (23) 
re X tana - Ve (t - T) cos A, 6h 


where t is the elapsed time since the detonation wave passed the apex of the 


cone, T = a = X/(u cos a) and the angle A = %+ 5, The slope of the contour 
D 

of the collapsing liner at any time t can be obtained from the Ir/5Z using 

equations (23) and (24) and some calculus. The time when a given element 

reaches the axis is found from equation (24) with r = 0 or 


X tan oa ; 


t-T. Vy cos A (25) 


and dr/ 2% evaluated at r = 0 is the tan g. Using the Taylor relationship, 
equation (14), to simplify we have finally 


sina +2 sin 6cosA-X sina (1-t : 
dane ae eee) 
cos a - 2 sin 6 sin A+ X sin a (tan A + tan 6) V/V 
oO 


Since 26 = gr - q and 2A = 8 + qg, equation (26) may be further 
simplified to yield, 


; + j 
sin B - X sina (1 -tanAt 
tan § = — + : — ee ; at 
cos ®8 + X sin a (tan A + tan 6) v/v. 
° 


The term sR represents the partial derivative of YG with respect to X. 


t 
The angle 8 > 8* since V a <0, i.e., the collapse velocity decreases 


from apex to base, and for cone angles (2q) which are not extremely large. 
Additional mathematical details and a more complete discussion can be found in 
Reference [2] and Eichelberger's Doctoral dissertation [4]. 


The reader should pay careful attention to the quadrant in which the 
angles calculated in equations (26) and (27) may fall. For angles in the 
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Figure 6. The Coordinate Directions. 
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second quadrant, minor modifications of the equation may be necessary and 
numerical algorithms that calculate trigonometric functions should be checked 
for the proper quadrant. 


The PER theory provides expressions for My mes Vv 6 and —, although 


3’ Ver 
the number of unknowns exceed the number of independent equations. 
Experimental evidence that verifies the model and the assumptions of the PER 
theory was given by Eichelberger and Pugh [5]. Also, Eichelberger [6] noted 

a systematic discrepency between the theoretical and experimental results. It 
was noted that the liner acceleration to the axis was assumed instantaneous in 
the PER model but the effect of a finite acceleration time would influence the 
value of g In 1962, Allison and Vitali [7], in another classic paper, used a 
radioactive tracer technique and obtained excellent agreement between the PER 
theory and experiment. 


Although Allison and Vitali obtained excellent agreement between the 
theory and experiment for a 42° unconfined copper cone, they were unable to 


get the PER theory to agree with the data from a 42° confined copper cone. 
Merendino and Jonas of the Ballistic Research Laboratory obtained similar data 


for a 60° unconfined copper cone which also failed to match the results of the 
PER model, G. H. Jonas recalls that Allison concluded that the PER model 
provides valid predictions for shallow angle, unconfined conical charges, but 
needs to be modified for wide angle or confined concial charges. 


The PER theory constitutes the basis of all one-dimensional shaped charge 
jet formation models. Many extensions have been made to the PER model and 
auxiliary equations have been added to facilitate and extend the solution. 


The plate bending angle or the Taylor angle (6) can be determined by the 
formula proposed by Richter [8] and Defournaux [9]: 


pek 


1 
=o + 
¢ 


where p and « are the density and thickness of the liner wall, respectively, 
and e is the explosive thickness driving the liner. K and +, are constants 


which are determined from the type of explosive used and the angle of 
incidence, i, which the detonation wave makes with the liner. Equation (28) 
is coupled with the PER theory in the one-dimensional jet formation codes BASC 
(70, 11] and DESC [12] and probably others. Defournaux [9] also presents a 
good discussion of the liner collapse and formation which is an excellent 
Supplement to the PER study. Kerdraon [13, 14] modified equation (28) by 
including terms to account for the effect of liner radius in imploding 
cylindrical liner charges. 


The PER theory assumes the liner element is instantaneously accelerated 
to the axis as depicted in Figure 7(a). The first order correction is the 
assumption that the liner acceleration is constant over a finite period of 
time, as first proposed by Eichelberger [6] and used in Reference [12]. The 
velocity then increases linearly over a short period, as shown in Figure 7(b), 
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to 


(a) Instantaneous Acceleration: 


(b) Constant Acceleration: 


to 


(c) Exponential Acceleration: 


Figure 7. Liner Acceleration Histories, Chou and Flis [15]. 
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until it reaches its final velocity or collapses on the axis. In Reference 
[12], the acceleration is given by 


pe? 
where Pad is the Chapman-Jouguet pressure of the explosive, po and ¢ are the 


liner density and thickness, respectively, and c is an empirical constant. 
Figure 7 was taken from Chou and Flis [15]. 


Possibly, a more realistic form of the liner velocity history is given as 
an exponential form by Randers-Pehrson [16], 


t-t 
V(t) = Vi [1 - exp{ - ——2 ]], 


as shown in Figure 7(c). This formula requires knowledge of the time 
constant, T, for which Chou, et al [17] recommend 


MV 
i] 


cJ 


where Mis the initial mass per unit area of the liner and cy and Cy are 
empirical constants. 

Considering the effects mentioned above, Chou, et al [17, 178] provide a 
more exact expression for the projection (Taylor)angle, 


v 
6a fo. 


1 1 
wi2t% *Ft VMS 
where the prime denotes differentiation along the liner. Randers-Pehrson [16] 
proposed a similiar formula. Both References [16, 17, 18] show better 


agreement with the hydrocode results than the Taylor angle (steady state) 
formula for 6, 


Models of this type, for acceleration to the liner axis during the 
collapse process, have been incorporated into the PER theory. These one 
dimensional computer programs use the analyses presented in this section to 
calculate the Taylor angle 5, the flow velocity V, the collapse velocity Vy? 


? 

the stagnation point velocity vy» the collapse velocity gradients V , the jet 
° 

velocity Vyr the slug velocity Vg? the collapse angle 8, the jet mass mye and 


the slug mass m+ The input values consist of the complete geometry and 
material density of the conical liner and the explosive charge. The mode of 
initiation of the charge and the type of HE (high explosive) must be known to 
determine the HE density and detonation velocity. The mode of initiation 


determines the shape of the detonation wave (i.e., plane or spherical) and the 
angle of incidence of the detonation wave to the liner wall. Also, the 
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relative distribution of the kinetic energy entering the jet and slug are 
given as 


dE. 2 
aa = cos* (a + § - 6/2) (29) 
dE 
and 
dE. 2 
ap = sin (6+ 6- /2) (30) 
respectively. 


The stretching jet radius is also calculated assuming the jet to be 
circular in cross-section. In addition, the material near the liner apex does 
not have sufficient time to reach its theoretical collapse velocity: and hence 
the earliest formed portions of the jet have lower velocities than the jet 
material formed behind them. This "piling up" of jet particles causes the 
formation of the inverse velocity gradient and this coalescence of particles 
is treated separately to calculate the jet tip velocity. All these effects 
are included in most jet formation codes such as BASC [10, 11] along with 
certain other parameters such as momentums, jet strains and strain rates. 


To calculate the strain, jet length and jet radius, it is necessary to 
relate the position of a jet element back to its original position in the 
liner. This was done by Carleone and Chou by using a simple Lagrangian 
coordinate definition [19]. The coordinate x is used for the axial liner 
position and the coordinate ¢ for the position of the jet as shown in Figure 
8. The jet position is given as 


5 (x, t) = 2 (x) + (t - t,) Vs (x), Gn 


where ty is the time when the liner element first arrives at the axis, Z (x) 
is the location of the formation, and v; (x) is the jet velocity. Note that 
vy is a funetion of a liner position but ¢ is a function of both liner 


position and time. Outside the inverse velocity gradient region, a liner 
element with a smaller x always proceeds (is ahead of, or has a larger ¢ than) 
an element with a larger x. The strain rate, 7, is given as 


4p - 8 (86/9x), 


ne SET Sanennticaiiby asa x-ree wate aE 
X > Xy Az (86/0x) (32) 


Equations (31) and (32), with the statements made earlier, allow the jet 
length and radius to be calculated. 

Special techniques are necessary to calculate the jet tip velocity as 
discussed earlier. The PER theory indicates that the jet velocity decreases 


IX-18 


LINER 


N 
\ — COLLAPSE PATH 
‘“ 
4+ +— 


x > 


2(x)}-———> 


—_ 


Figure 8. Relations Between the Liner Coordinate x, and the Jet Coordinate £ 
Chou and Flis [15]. 
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monotonically from tip to tail. Thus, the gradient of the jet velocity is 
negative with respect to the original liner position or positive with respect 
to the jet position. As mentioned earlier, liner elements near the apex 
collide on the axis before reaching their final collapse velocity, Vo This 


reduced collapse velocity results in a reduced jet speed, so that the jet 
velocity gradient with respect to the liner position is positive. This is 
termed the inverseevelocity gradient [20]. In the apex region, each jet 
element has a higher velocity than the one ahead of it, causing a "piling-up" 
of jet mass. This piled up mass forms the jet tip as verified experimentally 
[21]. For conventional conical charges, the first 30 to 40% of the liner from 
the theoretical apex forms the tip of the jet. The tip has a much larger 
radius than the rest of the jet. 


Figure 9 shows the inverse-velocity gradient with respect to X, the liner 
element position. The velocity of the jet tip is given as 


x, 
TIPy dm, dX 
J ax 


_ #0 
Veg * Voqp * ia , (33) 
dm. 
dx 
o a 


based on the conservation of momentum. 


Recently, Hirsch [22] extended the PER model to include non-steady 
effects at the collision point (or formation point) by assuming that relative 
to this point, the jet and slug velocities, v, and Vor are not equal to the 


ineoming liner velocity V as the PER theory assumed from the Bernoulli 
equation. 


Other studies involve modifications or extentions of the PER theory, such 
as Perez, et al [23] who modified PER to account for two-dimensional effects 
of the flow near the collapse axis, and Hirsch [24]. Leidel [25] provides an 
overview of the PER model and analyzes annular cutting charges, or cookie 
cutter warheads as first conceived by Glass, et al [26]. Harrison [27, 28] 
compares the simple one-dimensional models to hydrocodes and experimental 
data. These methods are also related to warhead design. 


Other shaped charge publications of interest include Hornemann and Senf 
[29] on wide angle cones, Van Thiel and Levatin [30], Trinks [31] for his 
shaped charge discussion and historical inferences, Rostoker [32], Simon and 
DiPersio [33], and the USSR studies by Novikov [34] and Titov [35]. These 
papers relate to jet formation and collapse and related topics pursuing the 
concepts we have reviewed so far. An excellent Soviet paper was published by 
Ivanova and Rozantseva [36] which is very long and very informative. It 
provides a detailed discussion and derivation of the PER theory and simple 
penetration models including experimental data and numerous flash radiographs 
or x-ray photographs of shaped charge jets. The USSR models are discussed in 
more detail later in the text. 
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Figure 9. A Typical Jet Velocity Distribution Curve Showing That the 
: Inverse-Velocity Gradient Region Near the Apex of the Liner 
Forms the Jet Tip Particle, Chou and Flis [15]. 
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Excellent flash radiography data and analysis thereof is given by 
Breidenbach [37] and by Heine-Geldern and Pugh [38]. Classical papers on jets 
and jet formation are given by Harlow and Pracht [39] and on collapsing 
cylinders by Koski, et al [40]. 


Zernow and Simon [41], investigated the behavior of jets from several 
different shaped charge materials. Aseltine [42] studied the effect of 
several different warhead asymmetries on the behavior and performance of 
shaped charge jets. Merendino, et al [43] altered a shaped charge liner, 
i.e., inhibited the collapse of the liner to obtain a massive hypervelocity 
pellet in lieu of a jet. Other studies regarding shaped-charge analyses are 
given in the recent papers by Robinson [44, 45] and Lee (46}. 


Kolsky, et al [47] studied the formation and penetration of conical 
shaped charge liners. Kolsky [48] continued his studies regarding the 
formation and penetration of hemispherical shaped charge liners. The 
formation and collapse of hemispherical liners were also studied via 
hydrocodes by Kiwan and Arbuckle [49], Arbuckle, et al [50], Chou, et al [51] 
Aseltine, et al [52], Chou, et al [53], and Lacetera and Walters [54]. 


Studies on the formation of hemispherical shaped charge liners were 
presented by Lee [46] Singh [55], Shepherd [56] and Grace [57]. A closed form 
analytical expression (analogous to the PER theory for conical liners) for the 
collapse and formation of hemispherical shaped-charge liners does not 
currently exist. 
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X. THE VISCO-PLASTIC THEORY 


Other studies involving the collapse and formation of shaped charge 
liners were advanced by Soviet researchers. Returning to the collapse of 
conical liners or wedges, the USSR used the results of Birkhoff, or perhaps 
Lavrenttev, as the basis for their studies, Lavrent'ev [1] having claimed 
prior discovery. In any case, the results of Birkhoff or Lavrent'ev were 
extended to include visco-plastic effects by Godunov, et al [2]. 


The visco-plastic model follows from References {1 - 7] and a detailed 
derivation and discussion of this model is given by Walters [8]. The basic 
viscoplastic theory [1 - 7] does not include transient effects. Typically, in 
the axisymmetric hydrocode models used at the BRL and other agencies, 
compressible flow is assumed and the material constitutive relationships are 
based on elastic-perfectly plastic, work hardening models or more 
sophisticated models, e.g. [9, 10]. The USSR analytical models typically 
assume an incompressible flow, but use a rate-dependent, visco-plastic 
material constitutive equation [2, 3, 11]. This equation is usually of the 


form 9 = - + é, where the stress 9%, is related to the strain rate é, by the 


yield stress ’ and the constant dynamic viscosity coefficient, v. 

Thus, the visco-plastic models require a knowledge of the dynamic 
viscosity coefficient and many USSR investigators have deduced viscosity 
coefficients from experimental measurements under shock loading conditions. 
Walters [11] summarizes many of the experimental viscosity values. Additional 
data related to viscosity measurements and the collapse of metallic bodies by 
high explosives are given by References [12 - 34]. Material viscosity values 
have also been experimentally determined in the US [11, 35], under shock 
loading conditions. The viscosity values deduced from various experimental 
measurements depend on many parameters, primarily strain rate, pressure and 
temperature. Thus, the viscosity coefficient is not constant. In fact, for 


solid metals, the dynamic viscosity may range from typically 10 to 10° Pa-sec 


(10° - 10° poise) depending on the strain rate, pressure and temperature and 
the experimental method used to measure }. Usually, the US measured viscosity 
values under shock loading conditions may be as much as two 
orders-of-magnitude lower than the USSR measured viscosity values [11 - 35]. 
The uncertainty in the US measured VU values is + a factor of ten. The USSR 
methods may even have a greater amount of uncertainty. Experimental 
determination of in the US, as performed by Chhabildas and Asay [35], 
involves measurement of the shock rise time via velocity interferometry. The 


maximum stress is taken to be proportional to the strain rate, SLE, and the 
maximum stress is determined from the Rayleigh line and the shock Hugoniot 
pressures. The strain rate is calculated as the strain induced by the shock 
divided by the shock rise time, which is resolution limited. The time 
resolution is limited to 1 to 3 nanoseconds (35). 


The USSR shaped-charge liner collapse model includes the influence of 
viscosity on the jet formation process. A jet formation criterion is derived 
which is based not on taking account of the compressibility effects or shock 


effects as in References [36] or [37], or on consideration of a critical Mach 
number as in Reference [38], but on taking into account the viscous properties 
of metals. 


Following Godunov, et al [2], an appropriate method is used to estimate 
the effect of viscosity in the plane problem (two-dimensional) of jet 
collisions. The fluid is incompressible, the motion irrotational and steady 
state, and the coefficient of viscosity 1s constant. Then the solutions of 
Euler equations automatically satisfy the Navier-Stokes equations of motion 
and the difference between the problems of ideal and viscous jet collisions is 
the conditions on the free surface of the jet. The flow fields are taken to 
be Newtonian. For the flow fields to agree in the ideal and viscous flow 
problems, some forces must be applied to the free surface in the viscous flow 
case [2]. Consideration of the influence of these forces on the flow 
resulting from the collision between jets of an ideal fluid will provide a 
measure of the influence of viscosity in the jet formation problem [2]. 


The formation of a reverse jet (slug) is possible only if the horizontal 
component of the viscous force, acting on the free surfaces of the reverse 
jet, is less that the force resulting from the symmetric collision of two 
plane fluid jets [2]. This inequality allows one to define a Reynolds number 
and a critical Reynolds number for jet formation. 


The visco-plastic shaped-charge jet collapse model derived by Godunov [2] 
relaxes to the classical US models reviewed earlier when the dynamic viscosity 
is zero. Basically, the visco-plastic model yields a lower jet velocity and a 
lower flow velocity than the US models. However, the visco-plastic model 
predicts a higher slug velocity than the US models. The USSR (visco-plastic) 
eriterion to form a coherent jet is that the Reynolds number, 


eee 

tusin 8 
Re = (1 _~ sing) me 
be greater than two. The wall thickness is denoted by t, v is the kinematic 
viscosity, 28 is the collision angle, ‘and u is the inviscid flow velocity as 
given by Defourneaux or the PER model for example. The Reynolds number 
greater than two criterion is based on experimental observations [2, 11]. 
Detailed USSR formulae for the jet strain rate, jet velocity, viscid flow 
velocity, and slug velocity, are given in References [2, 11]. Also, by 
setting Re = 2 in Equation 1, an expression for the critical flow velocity as 
a function of 8, t and y can be obtained. If the critical flow velocity is 
greater than the flow velocity, a coherent jet will not form. Simonov [6] and 
others [11 - 14, 18, 20, 21, 24, 27, 28, 31, 32] define collapse angles for 
the transition from the jet regime to the jetless regime for explosive bonding 
applications. These USSR investigators defined collapse angles below which 
only incoherent jets will form and collapse angles below which no jet at all 
will form. Walsh, et al [36] served as the basis for many of the USSR 
studies. 


For the transient flow (PER) case, the stagnation point velocity Wo) is 
given as 


sing 


Ve = D, [se - a)-sin(p ~ a - | ; (2) 


where Da is the speed of the detonation wave, 2a is the conical liner apex 


angle and 4 = 2618 the plate bending angle. Equation 2 is used in both the 
PER and visco-plastic models. The flow velocity is given by 


sing 


U=D, [eae +¢ =sins| (3) 


in the PER model. The visco-plastic flow velocity is 
U2 =U (1 - 2/Re)’# (4) 


from Reference [2] for Re > 2. Also, U. = U when the dynamic viscosity, y, is 


zero. 


2 


In the PER model, the jet velocity and slug velocity are, respectively, 


vy eV +U, (5) 
and 

Vs = Vo - U. (6) 

In the visco-plastic model, the jet and slug velocities are 

V5 = V, + Ups (7) 
and 

Vo 2 Vi, - U5. (8) 


In both models, 


, + v, s eV 


For steady state collapse, analogous formulae are available from Godunov 
and Birkhoff. 


e 


For both the steady state and transient version of the visco-plastic 
theory, the US theory can be obtained by setting the dynamic viscosity, u, or 
the kinematic viscosity, y = u/o equal to zero. Walters and Harrison [8] 


X-3 


discuss the PER and the visco=plastic models and comparison is made to the 
experimental data obtained by Allison and Vitali [39]. 


The strain rates calculated by the visco-plastic model (2] are of the 


order of 10°/s and show little variation with liner element position. Chou 
and Caleone [40] calculated strain rates for a 81 mm diameter precision charge 


with a copper conical liner with a 42° apex angle and a 1.9 mm wall thickness 
and found a larger variation of strain rate with a liner element position 
including liner elements in both compression and tension. The peak strain 


rate was 10°/s. Experimental measurements of the strain rate of collapsing 
liners are not available, but Bauer and Bless [41] measured the strain rate of 


exploding copper tubes to be about 10"/s and Walters deduced the strain rates 
of stretching jets to be of the order of 10'/s. . 


Comments on Jet Coherency 


The visco-plastic coherent jet or incoherent jet criterion is based on a 
critical Reynolds number being greater than 2 [2], or some other number [44]. 
The US criterion is inviscid, but compressibility effects are taken into 
consideration in the determination of the cohesiveness of the jet. 


Walsh, et al [36] studied the compressibility effect in the jet formation 
process as related to the oblique impact of explosively driven plates. The 
symmetric plate collision as viewed from a moving coordinate system, reduces 
to a flow field analogous to two impinging streams. Walsh, et al [36] 
concluded that jetting always occurs if the fluid is incompressible or if the 
collision velocity in the moving coordinate system is subsonic. For 
supersonic flow, jetting always occurs if B> Bor where B, is the critical 


turning angle for an attached oblique shock wave at the collision velocity. 
Walsh, et al [36] did not address jet cohesiveness per se, since their primary 
interest was explosive bonding. Cowan and Holtzmann [45] present additional 
flow criteria for explosive bonding. 


In the studies of Chou, et al [37, 46], jetting criteria for plane 
axisymmetric cases are presented along with a measure of the jet quality. 
From References [37] and [46], Chou states: 


1. For subsonic collisions (or the collision velocity V < C, the 
material bulk speed of sound) a solid coherent jet always forms; 

2. For supersonic collisions (V > C) jetting occurs if f& > Bys but the 
jet is not coherent. The angle 8, is the maximum angle that an attached shock 
wave can form at a prescribed supersonic velocity, V; and 


3. For supersonic collisions (V > Cc) but B< Bye a jet will not be 
formed. 


For shaped charge applications, the major criterion for jetting is that 
the formation process be subsonic. Otherwise, the jet will be incoherent and 
Spread out radially. 


The speed of sound referred to is usually taken to be the bulk speed of 
sound of the material as opposed to a longitudinal or transverse shear speed 
of sound. The appropriate speed of sound, as well as its exact value under 
the extreme pressures and temperatures encountered in the formation region is 
not well known. Also, the actual flow field is compressible and 
twoedimensional (at least) and hence the subsonic collision condition is 
useful only as a general principle. In practice, it has been observed that 
the value of the critical Mach number based on the static bulk speed of sound 
is about 1.2 for conical copper liners, Harrison [38, 47], and Walters and 
Harrison [8]. In other words, the liner velocity divided by the bulk speed of 
sound, i.e., the Mach number, could be as high as 1.2 and still achieve a 
cohesive jet for a conical, copper liner. A flow velocity of 4.8 km/s (1.23 
times the bulk speed of sound of copper) is the calculated flow velocity above 


which 20° copper liners have been observed to have incoherent jet tips [38]. 
Figure 1 shows a flash radiography photograph of a copper, coherent and 
incoherent jet. 
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Figure 1. Radiographs of Jets From Two Typical Conical Charges. 


(a) 40° Copper-lined Charge, Subsonic Collision, 
Coherent Jet, and 


(b) 20° Copper-lined Charge, Supersonic Collision, 
Incoherent Jet. 
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XI. EXPLOSIVE WELDING, SHAPING AND FORMING 


As mentioned in the first section, we will briefly address the concepts 
of explosive welding, explosive bonding and explosive forming. 


The explosive working of metals and concepts related to explosive 
effects, standoff systems, contact systems, shock physics and material effects 
are presented in an excellent introductory book by Rinehart and Pearson [1]. 


An explosive forming procedure is shown in Figure 7. A die assembly 
holding the work piece is positioned in a tank of water. A vacumm is drawn 
below the work piece and an explosive charge is detonated in the water. The 
work piece then assumes the shape of the die (if all goes well). 


A method for bulging cylindrical tubing is shown in Figure 2. The 
explosive charge is sealed in plastic and positioned in the cylindrical work 
piece filled with water. The ends are sealed to add confinement and reduce 
water splash. The charge is detonated within the eylinder held by a split die 
assembly. Figure 2 shows the bulging cylinder produced. 


In explosive welding or bonding applications, two or more plates are 
driven together at high velocities and surface jetting occurs. This results 
in a bond between the two driven metals and can result in a wave interface of 
high strength between the welded materials. If the plate impact velocities 
are too high (relative to the material properties), the interface material 
will jet as a shaped-charge and the bonding will be poor. If the impact 
velocities are too low, no interaction will occur at the interface. The 
velocity must be such that plastic deformation occurs at the interface and 
local surface jetting occurs. The requirements for welding to occur are [2]: 


a. the existence of a jet at the interface, and 


b. an increase in pressure, associated with the rapid dissipation of 
kinetic energy, to a sufficient level for a sufficient time to achieve stable 
interatomic bonds. The pressure is determined by the impact velocity and the 
time available for bonding is determined by the velocity of the collision 
point. 


The final geometry of the explosive-metal system will vary depending on 
the geometry and materials of the parts to be welded. The nature of the 
interface between the welded parts will vary depending on the materials welded 
and the nature of the explosive-metal interaction system. 


Figure 3 shows a surface cladding test arrangement. The top plate, or 
flyer plate, only is driven by the explosive force in order to weld a thin 
sheet of material onto a heavy sheet. This technique is useful in Plating a 
good structural material with a thin cladding to protect it from a corrosive 
or hazardous environment. 


Explosive welding is similiar to explosive cladding, but two flyer plates 
are used to form an explosively welded final product. Figure 4 depicts an 
explosively welded sandwich or a three-layer weld. 

Other examples of the explosive working of metals as well as the details 
associated with these processes can be found in Reference [1]. 
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Blazynski [2] provides a recent collection of papers dealing with many 
aspects of the explosive working of metals. Pearson, in Chapter 17 of 
Reference [2], describes the explosive hardening of austenitic manganese 
steels subjected to severe impact and abrasion. These steels are used in 
railroad frogs, rock crusher jaws, grinding mills and similar devices. To 
explosively harden a steel, a thin layer of explosive is detonated in contact 
with the surface to be hardened, usually at grazing incidence, The 
propagation of the shock through the material increases the metal hardness 
(BHN) throughout the material. This also increases the yield and tensile 
strength of the metal, making it more resistant to impact wear and with less “ 
tendency to deform. Duvall [3] discusses shock wave propagation in metals. 


Pearson, in Reference [2], discusses explosive compaction where explosive 
loads are applied either directly or through a loading system to compact 
powders. Two methods involve the use of gun powder cartridges and the use of 
explosives to drive pistons into the powder compaction chamber. Explosive 
compaction can produce parts from powders that cannot be produced by 
conventional pressing. Improved material properties with high press densities 
ean be produced. Explosive compaction is extremely valuable in the powder 
metallurgy field. 


Pearson, Chapter 1 and El-Sobky, Chapter 6 of Reference [2], discuss 
explosive welding. We have briefly touched on this area and we will conclude 
our discussion with some of the major points of Reference [2]. 


Explosive welding is a solid phase welding process in which high 
explosives are used to join the weld surfaces in a high velocity collision. 
This collision produces severe, localized plastic flow at the interface 
between the two surfaces. Explosive welding is used to weld metal 
combinations many of which cannot be welded by conventional means. 


The weld is usually formed by oblique impact of the plate surfaces with 
the weld, progressing from the apex along the collision interface by a process 
similar to the collapse of a shaped-charge liner. The parameters involved in 
explosive welding include the physical and mechanical properties of the metals 
to be welded, the type and amount of explosive used, the mode of initiation, 
the initial geometry of the weld operation and the type and geometry of buffer 
' sheets (if used) between the metal and the explosive. These parameters 
influence the collision angle, the impact velocity and the collision point 
velocity. 


It is generally agreed that jet formation at the collision point is 
essential for welding, e.g., El-Sobky, Chapter 6 of Reference [2]. This is 
because, at least in part, jetting produces chemically clean surfaces which 
allow interatomic bonding under the high pressure conditions. 


The jet-no-jet criterion, or the limiting condition for jet formation, 
discussed earlier, is based on the studies of Cowan and Holtzman [4], Walsh, 
et ai [5], Chou, et al [6] and Harlow and Pracht [7]. These studies relate to 
shock formation and critical collapse angles and are also applied to 
shaped-charge jets. The mathematical description for the explosive welding of 
oblique plates involves relationships similar to those used in jet formation 
theory and in the Gurney type equations. In addition to Gurney, impact 
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velocity equations are given by Deribas, et al [8] and others as reviewed by 
El-Sobky, Chapter 6 of Reference [2]. El-Sobky reviews the jet formation 
criterion as well as minimum impact pressures and standoff distance 
requirements which must be met. El-Sobky notes that the collision velocity 
and plate impact velocity should be less than the bulk speed of sound of 
either welding component. However, Wylie, et al [9] suggests that the bulk 
speed of sound may be exceeded by as much as 25% with satisfactory welds. It 
is interesting that this coincides with the jet-no-jet criterion observed by 
Harrison, et al [10,11,12] for copper jets, i.e., Mach number 1.23. Also, the 
visco-plastic critical Reynolds number criterion, discussed earlier, was 
derived from explosive welding experiments [13,74]. 


A major area of study in explosive welding is the interfacial periodic 
deformation or the interfacial waves. Often a regular wavy interface pattern 
appears, but the nature of the interfacial waves vary depending on the 
materials and system geometry used. 


When regular waves are produced, they are very similar to the surface 
deformations found on solid surfaces subjected to erosive action of high speed 
liquid films, jet drags, e.g., in pump and steam turbine blades, and at the 
air-sea interface, e.g., ocean waves. A review of the mechanisms for wave 
formation are given by El-Sobky in Chapter 6 of Reference [2], 


A good general discussion of explosive welding is given in Reference 
[15]. Jetting phenomena, the weld interface, and the explosive load and 
standoff distance required for various weld combinations are discussed. Other 
excellent sources of information include the USSR sources, given in the 
discussion on the viscoeplastic jet formation model and viscosity measurements 
(13,14]. Schroeder [16] recently edited a text on explosive welding and 
related areas. Deribas [17] authored a classic text on explosive welding and 
work-hardening. Duvall [3] presents an excellent supplementary discussion on 
shock propagation and transmission through solids. Sewell [18] also discusses 
explosive welding regarding jet-no-jet criterion. : 


This concludes the discussion of explosive working of metals. The 
information provided in this section is not complete, nor was it intended to 
be. The intent was merely to provide an overview of explosive welding and 
related fields. Our goal was to provide a brief introduction and a few 
fundamental references into this field and to note the similarity to shaped- 
charge jetting phenomenon. 


WATER LEVEL 


I EXPLOSIVE 
[] CHARGE 


HOLD DOWN WORK PIECE 


RING 


Figure 1. Explosive Forming, 


XI-4 


TOP DIE 


EXPLOSIVE 
CHARGE 


BOTTOM DIE 


BEFORE 


AFTER 


Figure 2. Bulging Cylinder Tube. 
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Figure 3. Surface Cladding: 
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Figure 4. Explosive Welding, a Three Layer Weld (Sandwich). 
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